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INTERACTIVE  GRAPHICS  PLOTTING  SYSTEM  (IGPS) 


1 .  INTRODUCTION 

IGPS  (Interactive  Graphics  Plotting  System)  is  a  computer- 
graphics  program,  written  in  BASIC  for  the  Tektronix  4051  desktop 
computer,  which  generates  high-quality  plots  of  two-dimensional 
data;  the  plots  can  be  used  for  slide  presentations  or  can  be 
inserted  into  reports. 

The  program  is  interactive  and  prompts  the  user  for  all 
the  information  required  to  construct  a  plot?  furthermore,  the 
program  is  structured  to  permit  rapid  changes  to  an  existing  plot. 
Therefore,  the  user  can  keep  modifying  the  plot  he  is  drawing  and 
can  preview  it  on  the  screen  of  the  Tektronix  4051  until  he  is 
satisfied  with  its  appearance.  Then,  simply  by  depressing  one 
key  on  the  keyboard,  the  user  can  exactly  reproduce  his  plot  on  a 
Tektronix  4662  x-y  plotter. 

Figure  1.1  is  a  sample  plot  produced  in  this  fashion  by 
IGPS  on  a  4662.  The  plot  was  drawn  directly  onto  a  piece  of 
report  paper  8-1/2  inches  high  by  11  inches  long  placed  in  the 
lower  left-hand  corner  of  the  4662.  The  figure  number,  title, 
and  page  number  were  then  added  with  a  conventional  typewriter . 

The  sample,  which  could  have  been  used  directly  to  make 
a  slide,  displays  some  of  the  features  of  IGPS;  an  unusual  axis 
system  (in  fact,  Rayleigh  Paper),  neat  tick  mark  and  axis  labeling, 
multiple  curve  marker  types,  legends,  and  curve  fitting.  Table 
1-1  is  a  summary  of  the  features  currently  available  with  IGPS. 
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Figure  1.1  Sample  Plot  Produced  by  IGPS 
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CUMULATIVE  FREQUENCY  CJO 


TABLE  1-1  SUMMARY  OF  IGPS 


Tick  Mark  Spacing: 


•  uniform 

«  user-controlled 

Tick  Mark  Labeling: 

•  uniform  with  automatic  decimal-point  selection 

•  user-controlled 

Frame  Types: 

•  axis 

•  grid 

•  border 

•  suppressed 

Standard  Paper  Types: 

•  linear  in  x  and  y 

•  linear  in  x,  logarithmic  in  y 

•  logarithmic  in  x,  linear  in  y 

•  logarithmic  in  x,  logarithmic  in  y 

•  user-definable 


Probability  Paper  Types: 

•  uniform  distribution 

•  normal  distribution 

•  '  lognormal  distribution 

•  Rayleigh  distribution 

•  Weibull  distribution 

#<  user-defined  distribution 

Input/Output  Devices: 

•  Tektronix  4051  Internal  Tape  Drive 
»  Tektronix  4924  Auxiliary  Tape  Drive 

•  Tektronix  4907  File  Manager  (Disc) 

Line  Types: 

•  solid 

•  dashed 

•  dashed-dotted 

•  dotted 
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TABLE  1-1  SUMMARY  OF  IGPS  (continued) 


Curve  Markers: 

•  crosses 

•  x-shaped 

•  squares 

•  triangles 

•  upside-down  triangles 
Curve  Fits  (Linear,  Least-Squares): 

•  estimate  slope  and  intercept 

•  estimate  slope,  intercept  fixed 

•  estimate  intercept,  slope  fixed 

Plotting  Devices: 

•  Tektronix  4051  Screen 

•  Tektronix  4662 

•  Tektronix  4662,  Option  31 


» 
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2. 


USING  IGPS 


Once  IGPS  has  been  loaded  into  the  memory  of  the  Tektronix 
4051,  it  needs  to  be  initialized.  This  can  be  done  by  typing  in 
RUN  followed  by  a  carrige  return.  IGPS  will  initialize  itself 
and  will  display  the  following  menu  on  the  screen  of  the  4051. 


LIMITS 

LABELS 

FORMAT 

10 

FITS 

! 

2 

3 

4 

5 

Ubt.K  l/Lr  IInADLL 

4051 

4662 

S&R 

FUNCTION 

MENU 

6 

7 

8 

9 

10 

Figure  2.1  Menu  for  IGPS 


The  10  rectangles  in  the  menu  correspond  to  the  10  user-definable 
keys  located  in  the  upper  left-hand  corner  of  the  keyboard  of  the 
4051.  These  10  keys  are  used  to  control  IGPS. 

If,  at  any  time  during  a  plotting  session  with  IGPS,  the 
user  wishes  to  examine  the  menu,  he  can  either  type  in  RUN  followed 
by  a  carriage  return  or  he  can  depress  user  key  number  10. 

Table  2-1  summarizes  the  action  of  the  10  user-definable 

keys . 


T 


t. 
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TABLE  2-1  SUMMARY  OF  USER  KEYS 


User  Key  Number  Function 

1  Establish  limits  over 

which  plotting  is  done. 

2  Inputs  axis  labels  and  title. 

3  Defines  axis  system. 

4  Inputs  data  curves  from 

peripheral  or  keyboard; 
saves  data  curves  on 
peripheral . 

5  Does  least-squares,  linear 

fits  on  data  curves. 
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Plots  graph  on  screen  of 
4051. 


7  Plots  graph  on  4662  x-y 

plotter. 

8  Stores  graph  on  periperal; 

retrieves  previously  saved 
graph  from  peripheral. 

8  Generates  x-y  data  from 

user-defined  parametric 
functions . 

10  Displays  menu  on  the 

screen . 


2.1  User  Key  Number  1;  Limits. 

This  key  is  used  to  establish  the  limits  between  which 
data  plotting  is  to  take  place  (in  other  words,  data  outside  the 
limits  will  not  be  plotted);  the  key  is  also  used  to  define  the 
distance  between  tick  marks  and  between  tick  mark  labels. 

After  the  user  depresses  key  number  1,  IGPS  will  prompt 

him  with 

UNIFORM  TICKS  ON  THE  X  AXIS  (1=YES,  2=NO): 

If  the  user  responds  with  a  1,  meaning  that  he  wants  uniformly 
spaced  tick  marks  and  labels,  IGPS  will  prompt  with 


« 
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XMIN : 

XMAX : 

MAJOR  TICKS: 

MINOR  TICKS: 

The  user  should  respond  to  each  prompt  with  a  suitable  number. 

The  response  to  XMIN  establishes  the  smallest  value  on  the  x  axis 
and  that  to  XMAX  the  largest.  The  response  to  MAJOR  TICKS 
establishes  the  distance  between  adjacent  tick  mark  labels  and 
the  response  to  MINOR  TICKS  the  distance  between  adjacent  tick 
marks.  The  sequence  to  generate  the  x  axis  of  Figure  1.1,  for 
instance,  is 

XMIN:  0 
XMAX:  96 
MAJOR  TICKS:  12 
MINOR  TICKS:  4 

In  no  case  should  XMIN  ever  be  the  same  as  XMAX  and,  for  best 
results,  MAJOR  TICKS  and  MINOR  TICKS  should  be  multiples  of  the 
difference  between  XMAX  and  XMIN;  furthermore,  MAJOR  TICKS  should 
be  a  multiple  of  MINOR  TICKS.  Results  are  unpredictable  if  these 
conditions  are  not  met.  In  fact,  IGPS  will  create  n  tick  marks, 
where 

n  =  1  +  [(XMAX  -  XMIN)  /  MINOR  TICKS]  , 

and  the  quantity  in  brackets  denotes  "the  largest  integer  not  to 
exceed,"  and  the  tick  marks  will  appear  at  points  x^  such  that 

Xj^  =  XMIN  +  MINOR  TICKS  *  (i-1),  i  =  l,...,n. 

Similarly,  IGPS  will  create  m  tick  mark  labels,  where 

m  =  1  +  [(XMAX  -  XMIN)  /  MAJOR  TICKS]  , 

and  the  labels  will  appear  at  points  1^  such  that 


1^  =  XMIN  +  MAJOR  TICKS  *  (i-1),  i=l,...,m. 

Clearly,  if  MAJOR  TICKS  is  the  same  as  MINOR  TICKS,  then  all  tick 
marks  will  be  labeled. 

Consider  the  sequence 

XMIN:  0 
XMAX:  1 

MAJOR  TICKS:  0.25 
MINOR  TICKS:  0.25 
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This  produces  the  following  x  axis 


0.00  0  *  25  0.50  0.75  1.00 


Figure  2.2  Sample  Axis  with  Uniform  Tick  Marks 
and  Labels 


Note  that  the  labels  are  drawn  centered  under  the  appropriate 
tick  mark,  and  that  all  labels  have  the  same  number  of  digits 
(two)  to  the  right  of  the  decimal  point.  Note  also  that  two 
digits  beyond  the  decimal  point  is  the  ideal  number  in  this 
situation,  in  the  sense  that  using  any  fewer  would  have  caused  a 
loss  of  accuracy  while  using  any  more  would  have  caused  unnecessary 
trailing  zeroes.  IGPS  calculates  and  uses  this  ideal  number 
automatically . 

If  the  user  does  not  like  uniformly  spaced  tick  marks 
and  tick  mark  labels  with  automatic  decimal-point  selection  on 
the  x  axis,  then  he  should  respond  with  a  2  to  the  prompt 

UNIFORM  TICKS  ON  THE  X  AXIS  (1  =  YES,  2=N0 )  : 

By  doing  so,  he  assumes  total  control  over  the  placement  and 
appearance  of  tick  marks  and  labels;  after  such  a  response,  IGPS 
prompts  with 

NUMBER  OF  TICK  MARKS  (AT  LEAST  2): 

LOCATION  OF  TICK  MARKS; 

NUMBER  OF  TICK  MARK  LABELS; 

LOCATION  OF  TICK  MARK  LABELS: 

The  meaning  of  these  prompts  is  obvious. 

At  least  two  tick  marks  are  required  since  the  first  one 
entered  is  assumed  to  be  the  minimum  value  along  the  x  axis  while 
the  last  one  is  assumed  to  be  the  maximum  value.  The  number  of 
tick  mark  labels,  however,  is  unrestricted  and  may  even  be  zero; 
if  it  is  zero,  the  fourth  prompt  is  not  displayed  since  it  is  not 
needed . 


The  locations  of  tick  mark  labels  must  all  be  entered  on 
one  line  and  must  be  separated  by  a  single  blank  or  by  a  comma. 
The  labels  will  be  plotted  by  IGPS  exactly  as  they  were  entered 
(in  other  words,  IGPS  will  not  calculate  the  number  of  places 
beyond  the  decimal  point  for  the  labels  but  will  use  as  many  as 
the  user  himself  typed  in). 
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For  example,  the  sequence 


NUMBER  OF  TICK  MARKS  (AT  LEAST  2):  19 

LOCATION  OF  TICK  MARKS:  .1  .2  .3  .4  .5  .6  .7  .8  .9 

123456789 

10 

NUMBER  OF  TICK  MARK  LABELS:  3 

LOCATION  OF  TICK  MARK  LABELS:  .1  1  10 

if  used  in  conjunction  with  a  logarithmic  x  axis,  will  produce 


l _ i _ i _ i  i  i _ i.  i-t-j - 1 _ i _ i — i — i — i— 

,  I  1  10 

Figure  2.3  Sample  Axis  with  User-Controlled  Tick  Marks 
and  Labels 


Once  the  user  has  supplied  the  specifications  for  the  x 
axis,  IGPS  will  prompt  for  the  y  axis  in  exactly  the  same  way  as 
it  did  for  the  x  axis. 

2.2  User  Key  Number  2:  Labels. 

This  key  is  used  to  establish  the  axis  labels  and  the 
title  of  the  plot.  If  the  user  does  not  want  these  items,  he 
should  respond  with  a  carriage  return  when  prompted  for  them. 

The  x-axis  label  is  plotted  below  the  plotting  area  and 
is  centered  between  the  x-axis  limits.  The  title  of  the  plot  is 
similarly  centered  but  is  plotted  above  the  plotting  area  and  is 
underlined  twice  (if  the  plot  has  no  title,  no  underlining  takes 
place).  The  y-axis  label  is  plotted  to  the  left  of  the  plotting 
area  and  is  centered  between  the  y-axis  limits.  When  plotted  on 
the  4662  x-y  plotter,  the  y-axis  label  will  be  written  along  a 
baseline  rotated  90  degrees  from  the  horizontal  (see  Figure  1.1); 
on  the  screen  of  the  4051,  however,  the  individual  characters  are 
plotted  on  the  horizontal,  with  the  first  character  above  the 
second,  the  second  above  the  third,  and  so  on  (as  in  Chinese  or 
Japanese )--this  is  done  because  the  character-generator  of  the 
4051  cannot  produce  rotated  characters.  Close  examination  of 
Figure  1.1  will  reveal  that  different  character  sizes  are  used  in 
various  parts  of  the  plot;  this  feature  is  also  only  available  on 
the  4662  because  the  character-generator  of  the  4662  can  produce 
characters  of  arbitrary  size,  unlike  the  generator  of  the  4051. 

2.3  User  Key  Number  3:  Format. 

This  key  is  used  to  control  the  basic  appearance  of  the 
graph:  it  determines  the  type  of  frame  which  will  be  used,  the 

speed  at  which  the  graph  will  be  drawn,  the  color  of  the  frame 
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and  labels  and,  more  importantly,  the  type  of  paper  which  will  be 
used . 


After  the  user  depresses  user  key  number  3,  IGPS  will 
ask  for  the  type  of  frame  the  user  wants.  There  are  four  possible 
answers:  zero,  meaning  no  frame  or  labels  are  to  be  drawn;  one, 

meaning  draw  an  axis  system;  two,  meaning  draw  a  grid;  and  three, 
meaning  draw  borders. 

The  first  option  is  useful  when  a  set  of  curves  is  to  be 
added  to  a  graph  because  it  avoids  re-plotting  most  of  the  graph. 

The  second  option  produces  axes  running  along  the  left- 
hand  side  and  the  bottom  of  the  plotting  area. 


j. 


■  L  L.  l-l 


-1 — I — I — I 


Figure  2.4  Sample  Axis  System 


14 


The  third  option  produces  a  grid  over  the  plotting  area 
each  tick  mark  is  indicated  with  a  solid  line. 


Figure  2.5  Sample  Grid 


The  fouth  option  produces  axes  and  tick  marks  along  the 
edges  of  the  plotting  area. 


t - 1 - 1 — i — i — i — m - 1 - 1 - 1 — i — i — i  T~r 


J _ I  l _ I _ I _ L  -I  JL_J _ I _ 1 _ I _ I _ I _ I _ l_L 


Figure  2.5  Sample  Borders 
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After  he  selects  a  frame  type,  the  user  will  be  prompted 
for  the  color  of  the  frame  and  labels.  He  must  reply  with  an 
integer  between  one  and  eight,  corresponding  to  the  eight  pens  in 
the  turret  of  the  4662  x-y  plotter  equipped  with  Option  31.  The 
choice  of  color  is  ignored  if  plotting  is  either  on  the  screen  of 
the  4051  or  on  a  standard  4662. 

The  user  is  then  asked  for  a  pen  speed,  to  which  he  must 
respond  with  an  integer  between  10  and  570.  The  answer  is  inter¬ 
preted  in  millimeters  per  second  and  is  the  maximum  speed  at 
which  the  4662,  Option  31,  plotter  will  draw  (the  response  is 
ignored  on  the  4051  or  the  standard  4662).  For  previewing,  the 
user  should  select  570  mm/s.  However,  for  best  results,  he  should 
select  10  mm/s  since  the  lower  speed  leads  to  much  sharper 
definition  of  vectors. 

The  final  prompt  associated  with  user  key  number  3  is 
for  the  type  of  paper  the  user  wants.  IGPS  will  display 


PAPER 

TYPES: 

X-AXIS 

Y-AXIS 

1. 

LINEAR: 

X 

Y 

2. 

XLOG : 

LOG ( X ) 

Y 

3. 

YLOG : 

X 

LOG ( Y ) 

4. 

LOG LOG: 

LOG ( X ) 

LOG ( Y ) 

5. 

USER: 

? 

p 

6  . 

UNIFORM: 

X/100 

Y 

7. 

NORMAL: 

N ( X/100 ) 

Y 

8. 

LOGNORMAL : 

N ( X/10U ) 

LOG ( Y ) 

9. 

RAYLEIGH: 

SQR(-2*LOG( l-X/100 ) ) 

Y 

.0. 

WEIBULL: 

LOG ( —LOG ( 1- X/100 ) ) 

LOG ( Y ) 

PAPER 

TYPE: 

to  which  the  user  should  respond  with  an  integer  between  one  and 
10,  corresponding  to  his  choice  of  paper.  As  indicated  by  the 
prompt,  the  user's  choice  is  in  fact  a  selection  of  two 
transformations,  or  functions,  one  function  for  the  x  axis  and 
one  for  the  y  axis.  All  data,  including  tick  mark  locations 
(but  not  tick  mark  labels)  related  to  the  x  coordinate  are 
transformed  via  the  x-axis  function  before  being  plotted,  as 
are  all  data  related  to  the  y  coordinate  via  the  y-axis  function. 
Since  the  reader  may  not  be  fimiliar  with  all  available  paper 
types,  we  will  describe  each  one,  with  its  associated  functions, 
separately. 
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2.3a  Linear  Paper.  This  is  the  most  commonly  used  type  of 
paper  and  is  in  fact  plain  graph  paper.  The  two  functions 
associated  with  this  paper  are  the  identity  functions  and  therefore 
have  no  effect. 


Figure  2.6  Example  of  the  use  of  Linear  Paper 


2.3b  Xlog  Paper.  This  paper  has  a  logarithmic  x  axis  and  a 
linear  y  axis;  therefore,  all  x-axis  data  are  transformed  via  the 
function 


Fa ( x )  =  In  x 

before  being  plotted.  Because  of  this  transformation,  there  are 
implications  for  the  meaning  of  the  slope  and  intercept  of  any 
straight  line  which  the  user  might  choose  to  fit  to  his  data 
while  using  xlog  paper  (Section  2.5,  User  Key  Number  5).  In 
particular,  if  we  let  a  be  the  slope  of  such  a  line  and  if  we  let 
b  be  its  intercept,  and  if  the  line  appears  to  fit  the  user's 
data  plotted  on  xlog  paper,  then  an  appropriate  model  of  the 
user ' s  data  is: 

y  =  a  In  x  +  b  . 

Though  this  may  be  perfectly  obvious  in  this  case,  we  have  chosen 
to  emphasize  the  point  here  because  the  situation  becomes  much 
more  complex  with  other  types  of  paper  (transformations).  In 
general,  if  we  denote  the  two  axis  transformations  by  Fa(x)  and 
Fb(y),  then  any  straight  line  used  with  those  transformations  is 
in  fact  a  model  of  the  form: 

Fb(y)  =  a  Fa ( x )  +  b  . 
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If  the  user  wants  the  model  in  simpler  form,  with  only  the  y  term 
appearing  to  the  left  of  the  equal  sign,  he  will  have  to  apply  the 
inverse  of  F5  to  both  sides  of  the  equation;  this  point  is 
illustrated  explicitly  in  the  next  section. 

2.3c  Ylog  Paper.  This  paper  has  a  linear  x  axis  and  a 
logarithmic  y  axis.  To  illustrate  the  use  of  this  paper,  consider 
the  data  in  Table  2-2. 


TABLE  2-2  VALUES  OF  AIR  DENSITY  AS  A  FUNCTION  OF  ALTITUDE 


Altitude  ( km) 

0 

1 

2 

3 

4 

5 

6 

Air  Density  (kg/m^) 

1.22 

1.11 

1.01 

0.91 

0.82 

0.74 

0.66 

Altitude  (km) 

7 

8 

9 

10 

11 

12 

13 

Air  Density  (kg/m^) 

0.59 

0.53 

0.47 

0.41 

0.36 

0.31 

0.27 

Altitude  (km) 

14 

15 

16 

17 

18 

19 

20 

Air  Density  (kg/rn3) 

0.23 

0.19 

0.17 

0.14 

0.12 

0.10 

0.09 

From  physical  reasoning,  or  for  some  other  reason,  we  decide  that 
air  density  is  an  exponential  function  of  altitude  (in  other 
words,  that  the  logarithm  of  air  density  is  a  linear  function  of 
altitude).  To  test  this  hypothesis,  we  plot  the  data  on  ylog 
paper  and  fit  a  straight  line  with  user  key  number  5. 
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ALTITUDE  CKILOMETERS^ 


Figure  2.7  Example  of  the  Use  of  Ylog  Paper 


The  plot  confirms  our  hypothesis.  The  slope  and  intercept  of  the 
line  calculated  by  IGPS  are: 

♦ 

a  =  -0.13, 

and 

b  =  0.34  . 

Therefore,  our  model  of  air  density  (p)  versus  altitude  (h)  is: 

In  p  =  -0.13  h  +  0.34  , 
or 

p  =  1.40 

Note  that  this  model  does  not  go  through  the  data  point  at  zero 
altitude.  Suppose  we  should  want  it  to,  in  other  words,  suppose 
we  should  like  a  model  of  the  form: 

p  =  1.22  e-kh  . 

what  is  the  best  value  of  k  (in  the  least-squares  sense)?  Taking 
logarithms  produces: 


In  p  =  -kh  +  0.20  . 
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To  determine  k,  we  depress  user  key  number  5  and  request  a  straight- 
line  fit  with  a  forced  intercept  of 

b  =  0.20  . 

IGPS  will  calculate  the  slope  of  that  line  as: 

a  =  -0.12 

It  is  clear,  therefore,  that  the  value  of  k  we  are  looking  for  is: 

k  =  0.12  , 

and  that  the  best  model  which  goes  through  the  value  of  air 
density  at  zero  altitude  is: 


P 


1.22 


e-0.12h 


2.3d  Loglog  Paper.  This  paper  has  a  logarithmic  x  axis  and  a 
logarithmic  y  axis.  Therefore,  straight-line  fits  to  data  plotted 
on  this  paper  are  models  of  the  form: 

In  y  =  a  In  x  +  b  . 

The  user  is  cautioned  that  whenever  he  is  using  a 
logarithmic  axis,  be  it  horizontal  or  vertical,  he  must  avoid 
zero  or  negative  values  along  that  axis,  both  in  the  data  being 
plotted  and  in  the  choice  of  tick  marks.  This  is  because  the 
logarithm  of  a  non-positive  number  is  undefined. 

2.3e  User-Definable  Paper.  This  is  one  of  the  most  powerful 
features  of  IGPS  because  it  allows  the  user  to  define  his  own 
transformations  for  the  axes  and  therefore  provides  an  unlimited 
supply  of  paper  types.  If  you  want  to  use  this  feature,  you  must 
define  your  transformations  before  you  depress  user  key  number  3 
and  you  must  then  select  user-definable  paper. 

To  define  your  transformations,  you  must  replace  lines 
4590  and  4600  in  the  program  with  your  own  functions.  The  trans¬ 
formation  for  the  x  axis  must  appear  on  line  4590  and  must  be  of 
the  form 

4590  DEF  FNA ( X ) = 

where  you  supply  the  information  to  the  right  of  the  equal  sign. 
Line  4600  is  for  the  y  axis  and  must  be  of  a  similar  form 

4600  DEF  FNB ( Y ) = 

To  illustrate  the  power  of  this  feature,  consider  the 
data  appearing  in  Table  2.3. 


20 


TABLE  2-3  SAMPLE  DATA 


X 

3.00 

2.70 

2.30 

1.80 

1.60 

1.27 

1.10 

y 

.225 

.241 

.267 

.304 

.321 

.351 

.365 

Plotting  the  data  on  ordinary,  linear,  graph  paper 
reveals  nothing  unusual  about  them. 


Incidentally,  the  source  of  the  data  and  the  units  associated  with 
them  do  not  concern  us  here.  What  does  concern  us  is  that  a  model 
of  the  form 

k 

Y  =  _ 

f- 

has  been  suggested  for  the  data.  To  test  this  model,  we  type  in 

4590  DEF  FNA{ X ) =-l/SQR(X ) 

4600  DEF  FNB( Y ) =Y 

we  select  user-definable  paper  and  we  ask  for  a  curve  fit  with  a 
forced  intercept  of  b=0. 
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0.5 


T 


0*1 

0.0 


2  3  4 

Figure  2.9  Sample  Data  Plotted  on  User-Defined  Paper 


IGPS  calculates  the  slope  of  the  line  as  a  =  -0.40.  This  means 
that  one  possible  model  of  the  data  is: 


y 


0.40 


Note  that  the  x  axis  on  user-definable  paper  starts  at  1,  whereas 
the  x  axis  on  ordinary  paper  starts  at  zero.  This  difference  is 
necessary  because  the  transformation  of  zero  on  user-definable 
paper  would  lead  to  an  undefined  number  (in  other  words, 


-1 

-  =  —DO  )  . 

formation  as 


Note  also  that  we  defined  the  x-axis  trans- 


Fa(x) 


FaU) 


1 

V*” 

i 


rather  than  as 


The  minus  sign  is  required  to  make  the  transformed  value  of  the 
minimum  along  the  x  axis  smaller  than  the  transformed  value  of 
the  maximum  (in  other  words,  if  the  minus  sign  had  not  been 
included,  we  would  have  had  Fa(l)  >  Fa(4),  which  is  not  allowed 
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on  the  Tektronix).  The  minus  sign  has  to  be  taken  into  account 
when  interpreting  the  results  of  the  curve  fitting  (which  was 
done  in  constructing  the  final  model  above). 

The  final  point  to  make  about  this  example  is  that  the 
value  of  the  free  parameter  in  our  model  could  have  been  read 
directly  off  the  graph  at  the  point  where  the  fitted  line  intersects 
the  vertical  line  x=l.  One  way  to  see  this  is  to  note  that,  in 
our  model,  y=k  when  x=l. 

What  this  example  has  illustrated  is  that  user-definable 
paper,  in  conjunction  with  curve-fitting,  can  be  used  to  test 
the  adequacy  of  a  large  class  of  models  in  representing  given 
data  (though  no  mention  of  this  has  yet  been  made,  it  will  be 
shown  in  Section  2.5  that  IGPS,  when  curve-fitting,  generates  a 
quantitative  esitmate  of  the  goodness-of-f i t  of  the  model;  this 
estimate  is  known  as  ) .  As  another  example,  consider  the  data 
in  Table  2-4. 


TABLE  2-4  MORE 

SAMPLE  DATA 

X 

l 

2 

3 

4 

y 

1.01 

2.28 

3.92 

6.13 

One  proposed  model  for  these  data  is: 

x 

y  =  -  • 

1.1  -  kx 

By  inverting  both  sides  of  the  equation  and  multiplying  by  -1 , 
this  model  can  be  transformed  into: 

1  1.1 

—  1,1  =  —  +  k  . 

y  x 

Therefore,  we  can  check  the  adequacy  of  the  original  model  by 
typing 


4590  DEF  FNA ( X ) =-l/X 
4600  DEF  FNB ( Y ) =-l/Y 


and  then  selecting  user-defined  paper  and  doing  a  curve  fit  with 
a  forced  slope  of  a  =  1.1. 
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Figure  2.10  Another  Example  of  User- Defined  Paper 


The  plot  shows  that  the  model  is  a  good  one  and  as  IGPS 
calculates  the  intercept  of  the  line  as  b  =  0.11,  we  know  the 
value  of  the  free  parameter  k.  Therefore  our  model  becomes: 


x 

Y  =  - 

1.1  -  O.llx 

This  example  illustrates  that  a  particular  model  may  need  to  be 
manipulated  before  it  can  be  used  to  define  a  paper  type.  in 
general,  to  define  a  paper  type,  a  model  will  have  to  be  put  into 
the  form: 


Fb(y)  =  a  Fa(x)  +  b  , 

where  Fa  and  F^  are  functions  of  x  and  y  and  where  a  and  b  are 
functions  of  the  free  parameters  of  the  original  model. 

The  paper  types  described  so  far  have  been  suitable  for 
sets  of  two-dimensional ,  x-versus-y,  data.  we  now  describe 
several  types  of  paper  suitable  for  testing  whether  a  set  of  data 
is  from  a  particular  probability  distribution. 


24 


2.3f  Uniform  Probability  Paper.  Let  us  imagine  that  we  have 
a  sample  of  n  points,  d £ ,  i  =  l,...,n,  from  a  uniform  distribution 
on  the  interval  (A,B).  Suppose  further  that  the  cumulative 
fraction  of  the  data  associated  with  the  itft  point  is  q^, 
i=l,...,n.  The  q^'s  can  be  determined  in  two  steps:  (1)  sort 

the  di'c  into  increasing  order  and  (2)  calculate  q^  as: 

i 

q  i  =  “  "  f  i=  1 ,  .  . .  ,  n  . 

1  +  n 


Since  the  d^'s  are  uniformly  distributed,  it  must  be  the  case 
that 

qi  =  F(di) ,  i=l». . . »n  , 

where  F  is  the  cumulative  distribution  function  (c.d.f.)  of  the 
uniform  distribution  on  (A,B): 


It  follows,  therefore,  that 

di  =  F-1(qi)  ,  i=l , 


x  <  A, 

A  x  S.  B,  and 
x  >  B. 


n 


r 


where  F~1  represents  the  inverse  of  F  (the  existence  of  F~1  can 
be  guaranteed  by  restricting  the  domain  of  F  to  the  interval 
A  _£  x  <.  B)  .  In  other  words, 


d^  (  B— A)  qi  "t  A,  1-1 ,  .  .  .  ,  n  , 

so  that  if  we  plot  the  d^'s  against  the  q^'s,  we  should  be  able 
to  fit  a  straight  line  to  the  plotted  points;  what  is  more,  if 
the  slope  of  the  line  is  a  and  its  intercept  is  b,  we  can  calculate 
the  parameters  of  the  distribution  from  a  =  B  -  A  and  b  =  A;  these 
two  equations  in  fact  imply  A  =  b  and  B  =  a  +  b. 

Consider  the  data  in  Table  2-5. 
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TABLE  2-5  UNIFORM  NUMBERS 


5.39 

5.87 

6.41 

5.00 

6.76 

6.01 

6.62 

5.64 

5.41 

6.28 

6.87 

5.65 

5.81 

6.50 

5.17 

6.41 

6.83 

6.74 

5.22 

5.16 

To  check  the  uniformity  of  the  data,  we  sort  them  in  increasing 
order,  generate  the  associated  cumulative  fractions,  and  plot 
them  after  asking  for  a  curve  fit  ( IGPS  will,  if  requested,  do 
the  sorting  and  the  generating  of  cumulative  fractions;  note 
that  the  fractions  are  expressed  as  percentages,  however).  The 
plot  provides  convincing  evidence  of  the  uniformity  of  the  num¬ 
bers;  furthermore,  the  curve  fit  produces  a  slope  a  =  2.24  and 
an  intercept  b  =  4.87.  We  conclude  that  the  numbers  are  from  a 
uniform  distribution  on  (4.87,  7.11). 


CUMULATIVE  FREQUENCY  CJO 


Figure  2.11  Example  of  the  Use  of  Uniform  Probability  Paper 
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2.3g  Normal  Probability  Paper.  This  paper  can  be  used  to 
check  if  a  given  sample  of  data  is  from  a  normal,  or  gaussian, 
distribution . 

Imagine  that  we  have  a  sample  of  n  points  d-^ ,  i=l,...,n, 
from  a  normal  distribution  with  unknown  parameters  m  and  s. 

We  can  transform  the  d^ 1 s  to  standard  form  by  setting 

d^  -  m 

z ^  ,  l-l, ...  ,n  . 

s 


Furthermore,  since  the  d^'s  are  normally  distributed,  then  the 
zi's  are  distributed  according  to  the  standard  normal  distribution, 
and  it  must  therefore  be  the  case  that 


q  i  F  (  z  i  )  ,  i  1  ,  .  .  .  ,  n  , 


where 
the  c 


the 

d.f . 

F  (  X ) 


q^'s  are  the  cumulative  fractions  of 
of  the  standard  normal  distribution: 
x 


e""-*-/^  dt. 


the  data  and  F 


is 


It  follows  that 

z^  =  F  (q^)  ,  i=l  ,  .  .  .  ,  n  , 
and  from  the  definition  of  z-^  we  have 


d  ^  -  m 


s 


( q i ) ,  i=l 


n 


d^  =  s  F  1  (q^)  +  in,  i=l,...,n  . 

Therefore,  if  we  plot  d,  against  F-1(qi),  i=l,...,n,  the  result  will 
approximate  a  straight  line  whose  slope  and  intercept  can  be  used 
directly  to  estimate  m  and  s.  Consider  the  data  in  Table  2-6. 


TABLE 

2-6 

NORMALLY 

DISTRIBUTED 

DATA 

( FREQUENCY) 

Frequency 

.12 

.22 

.29 

00 
r— 1 

• 

.15  .02 

Value 

3.0 

3.5 

4.0 

4.5 

5.0  5.5 
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The  same  data,  in  cumulative  form,  appears  in  Table  2-7. 


TABLE 

2-7 

NORMALLY 

DISTRIBUTED 

DATA  (CUMULATIVE) 

Cumulative 
Frequency  (%) 

12 

34 

63 

81 

96 

98 

Value 

3.0 

3.5 

4.0 

4.5 

5.0 

5.5 

We  enter  the  data  in  Table  2-7  into  IGPS  after  selecting 
normal  paper  and  we  do  a  curve  fit. 


CUMULATIVE  FREQUENCY  OO 


Figure  2.12  Example  of  the  Use  of  Normal  Probability  Paper 


Estimates  of  the  mean  and  standard  deviation  provided  by  the 
curve  fitting  are  m  =  3.82  and  s  =  0.75.  These  two  quantities 
can  also  be  estimated  from  the  graph:  the  mean  is  the  point  at 
which  the  line  intersects  the  50*-h  percentile  while  the  standard 
deviation  is  the  difference  between  this  value  and  the  value 
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corresponding  to  a  cumulative  fraction  of  0.8413.  These  assertions 
can  be  verified  as  follows:  our  model  is  of  the  form 


d  =  sF  3  ( q)  +  m  , 

and  therefore  d  =  m  where  F“l(q)  =  0  and  F“l(q)  =  0  when  q  = 

0.50.  To  find  s,  note  that  s  =  d  -  m  when  F-3(q)  =  l  and  this  is 
true  when  q  =  0.8413.  Another  way  to  see  this  is  to  note  that 
the  mean  and  median  of  the  standard  normal  distribution  are  equal 
while  one  standard  deviation  about  the  mean  spans  68.26  percent 
of  the  distribution  (and  50  percent  plus  one-half  of  68.26  per¬ 
cent  equals  84.13  percent). 


Note  that  the  x-axis  tick  mark  labels  on  the  graph  avoid 
the  values  zero  and  100.  This  should  always  be  done  when  using 
normal  paper  because  F“1(0)  =  -°o  and  F”l(l)  =  +°°  . 

Finally,  a  remark  about  F-3  itself:  since  F“1  cannot  be 
expressed  in  closed  form  (that  is  to  say  in  terms  of  sums  and 
products  of  functions  available  on  the  Tektronix  4051),  IGPS  uses 
an  approximation  of  F“1  due  to  Hastings  (reference  1,  page  192). 


The  approximation  is: 

F_1(q)  = 


where 


jG(l-q)  , 
(  — G  (  q  )  , 

G  (  q )  =  z 


z  = 


In 


q  >  *5  , 
q  1  •  5  , 

h(  z ) 
s(  z ) 

1 

3 


and 


h(z)  =  2.515517  +  0.802853  z  +  0.010328  z2  ,  and 
s ( z )  =  1  +  1.432788  z  +  0.189269  z2  +  0.001308  z3. 


It  is  claimed  that  the  error  made  in  using  the  approximation  is  never 
greater  than  0.0004  standard  deviations. 


2.3h  Lognormal  Probability  Paper.  Lognormal  paper  is  similar 
to  normal  paper  in  all  respects  save  one:  the  y  axis  of  lognormal 
paper  is  logrithmic  whereas  the  y  axis  of  normal  paper  is  linear. 
Therefore,  lognormal  paper  can  be  used  to  test  whether  a  sample 
of  data  is  from  a  lognormal  distribution,  that  is  to  say,  that 
the  logarithm  of  the  sample  is  normally  distributed. 

2.3i  Rayleigh  Probability  Paper.  The  Rayleigh  distribution 
has  c .d . f . 
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q  =  F(d)  =  1  -  e  2  \s ) 
This  can  be  readily  transformed  into 


d  =  s  y-2  In  (1  -  q)  . 

Therefore,  given  a  sample  of  n  points  d^,  i=l,...,n,  with  cumulative 
fractions  qj_,  i=l,...,n  we  can  test  whether  the  sample  is  from  a 
Rayleigh  distribution  by  plotting  d^  against  yj -2  Tn  ( 1  -  q^ )  , 

i=l,...,n.  Furthermore,  we  can  estimate  s  by  fitting  a  line  with 
intercept  b  =  0  and  slope  a  to  the  transformed  data  and  setting 
s  =  a.  The  data  in  Table  2-8  will  illustrate  the  procedure. 


TABLE  2-8  NUMBERS  FROM  A  RAYLEIGH  DISTRIBUTION 


320 

470 

16 

619 

356 

547 

264 

204 

429 

703 

265 

305 

501 

128 

468 

662 

608 

143 

121 

524 

We  input  the  data  to  IGPS ,  ask  IGPS  to  sort  it  and  generate 
the  cumulative  fractions,  ask  IGPS  for  a  curve  fit  with  a  forced 
intercept  of  b=0,  and  plot  the  results  on  Rayleigh  paper. 


CUMULATIVE  FREQUENCY  CO 


Figure  2.13  Example  of  the  use  of  Raleigh  Probability  Paper 
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One  estimate  of  the  parameter  s  is  provided  by  the  curve 
fit  as  s  =  314.  Another  is  provided  directly  from  the  graph  at 
the  point  where  the  line  intersects  the  3 percentile.  This  is 
the  case  because  our  model  is: 


d  =  s  yj  -  2  ln(l-q)  , 

and  therefore  d  =  s  when  -J -2  In  ( 1-q)  =  1,  which  is  the  case 
when  q  =  39%. 

Note  that  the  value  100  is  not  used  as  an  x-axis  tick 
mark  label;  this  is  because  when  q  =  1,  1^  (1-q)  =  +°°  • 

2.3j  Weibull  Probability  Paper.  The  Weibull  distribution 
has  c.d.f. 

c 

q  =  F(d )  =  1  -  e“kd,  d  >  0,  k  >  0. 

This  can  be  transformed,  by  taking  logarithms  twice,  into 

1  1 

In  d  =  -  In  [-In  (1  -  q)  ]  -  - -  In  k. 

c  c 

Weibull  paper  is  ruled  according  to  the  function 

Fa( x)  =  In  [-In  (1  -  x) ] 

along  the  x  axis  and  according  to  the  function 

Fb(y)  =  In  y 

along  the  y  axis.  Therefore,  given  a  set  of  data  d-[  ,  i  =  l,...,n, 
from  a  Weibull  distribution  with  parameters  c  and  k,  with  associated 
cumulative  fractions  q^ ,  i=l,...,n,  then,  if  we  plot  d^  against 
qi  ,  i  =  l,...,n,  on  Weibull  paper,  the  result  should  approximate  a 
straight  line.  Furthermore,  if  the  slope  of  the  line  is  a  and 
its  intercept  is  b,  then  estimates  of  c  and  k  are  given  by 

1 

c  =  -  ,  and 

a 

k  =  e~b/ a 

Note  that,  because  of  the  definition  of  Fa,  values  of  zero  and  100 
should  be  avoided  on  the  x-axis  (in  fact,  Fa(0)  =  -«  and  Fa(l)  =  +») 
similarly,  non-positive  values  should  not  be  used  along  the  y  axis. 

Before  leaving  this  section,  we  should  point  out  that 
the  user  can  define  his  own  probability  paper  by  using  the  ideas 
developed  in  Section  2.3e.  For  instance,  suppose  the  user  is 
interested  in  constructing  paper  for  the  exponential  distribution. 
Now  the  exponential  distribution  has  c.d.f. 
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q  =  F(d)  =  1  -  e_cd,  c  >  0. 


By  rearranging  and  taking  logarithms,  this  can  be  converted  to 

1 

d  =  -  ‘  '  In  (1  -  q)  . 

c 

Therefore,  to  construct  exponential  paper,  the  user  would  type  in: 
4590  DEF  FNA(X)=-LOG ( l-X/100 ) 

4600  DEF  FNB ( Y ) =Y 

and  would  then  select  user-definable  paper.  Note  that  the  argu¬ 
ment  of  the  function  FNA  is  divided  by  100;  this  is  done  to  allow 
the  x-axis  tick  marks  and  labels  to  be  expressed  as  percentages. 

We  leave  to  the  reader  the  problem  of  determining  what 
kind  of  curve  fit  to  use  with  exponential  paper  and  how  to 
estimate  the  parameter  c  from  the  results  of  the  curve  fit. 

The  overall  strategy  for  creating  distribution  paper, 
however,  is  clear:  the  c.d.f.  of  the  distribution  of  interest  is 

manipulated  into  the  form 

Fb(d)  =  a  Fa(q)  +  b, 

where  Fa  and  Fb  are  functions  of  the  cumulative  frequency  of  the 
data  and  of  the  data  itself  and  a  and  b  are  functions  of  the  free 
parameters  of  the  distribution. 

2.4  User  Key  Number  4:  Input  and  Output. 

This  key  is  used  either  to  input  data  curves  into  IGPS 
or  to  save  previously  entered  curves  from  IGPS  onto  a  Tektronix 
peripheral . 

2.4a  Input.  After  the  user  chooses  to  input  data  into 
IGPS,  he  will  be  prompted  with 


DATA  TYPE  (1  =  Y  ONLY,  2  =  X  &  Y)  : 

A  response  of  1  indicates  that  the  user  is  only  going  to  input  y- 
coordinate  data  whereas  a  response  of  2  indicates  that  he  is 
going  to  input  both  x-axis  and  y-axis  data.  A  response  of  1  is 
also  a  request  to  IGPS  to  do  two  things:  (1)  sort  the  y-axis 
data  in  increasing  order  and  (2)  generate  x-axis  data  as: 
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where  n  is  the  number  of  curves  in  the  data  set  and  m^  is  the 
number  of  points  in  the  ith  curve.  This  feature  is  useful  with 
probability  paper;  see  sections  2.3f  through  2.3j. 


The  user  is  then  prompted  for  the  device  he  wants  to 
read  the  data  curves  from.  There  are  four  possible  responses; 
one,  meaning  the  Tektronix  4907  File  Mangager,  a  disc  system; 
two,  the  Tektronix  4924  Auxiliary  Tape  Drive;  three,  the  Tektronix 
4051  Internal  Tape  Drive;  and  four,  meaning  the  user  wants  to 
input  data  into  IGPS  by  typing  it  on  the  keyboard  of  the  4051. 


The  user  is  then  prompted  for  the  number  of  curves  he 
wants  to  input  and  for  the  number  of  points  in  the  curve  having 
the  largest  number  of  points  (in  other  words,  for  the  number  mg 
where  mg  =  max  [  m^,  i  =  l,...,n  ]  ,  where  n  is  the  number  of 

curves  being  input  and  m^  is  the  number  of  points  in  the  i  ^ 
curve ) . 


If  the  user  selected  the  4907  as  the  device  to  use,  he 
is  then  prompted  for  the  names  of  the  files  on  which  his  data 
reside;  if  he  selected  either  the  4924  or  the  internal  tape  drive, 
he  is  then  prompted  for  the  numbers  of  the  files;  and  if  he 
selected  the  keyboard,  he  is  prompted  for  the  number  of  points 
and  for  the  values  of  the  data  points  in  each  curve  (furthermore, 
if  the  user  selected  y— only  data,  he  is  only  prompted  for  the  y— 
axis  values,  whereas,  if  he  selected  x  &  y  data,  he  is  prompted 
for  both  x-axis  and  y-axis  values). 

When  IGPS  has  read  in  all  the  data,  it  asks  for  the  type 
of  curve  marker  to  be  used  with  each  curve.  These  are  nine 
different  curve  markers  available: 
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Figure  2.14  Available  Curve  Markers 
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Interesting  effects  can  be  achieved  by  plotting  the  same 
curve  more  than  once,  with  a  different  curve  marker  each  time. 

IGPS  next  prompts  for  the  legend  the  user  wants  for  each 
curve;  if  the  user  does  not  want  a  legend,  he  should  respond  with 
a  carriage  return.  The  legend  for  a  particular  curve  is  plotted 
slightly  to  the  right  of  the  last  data  point  in  the  curve.  This 
feature  was  used  to  create  the  symbols  one  through  nine  in  the 
plot  above. 

Finally,  IGPS  prompts  for  the  color  the  user  wants  for 
each  curve.  Allowable  responses  are  integers  between  one  and 
eight;  these  refer  to  the  eight  pens  in  the  turret  of  the  Tektronix 
4662  x-y  plotter  equipped  with  Option  31.  The  responses  are 
ignored  if  plotting  is  either  on  the  screen  of  the  4051  or  on 
the  standard  4662. 

If  the  user  selected  y-only  data,  there  may  now  be  a 
slight  pause  as  IGPS  sorts  the  y-axis  data  and  generates  the  x- 
axis  data.  The  y-axis  data  are  sorted  in  place  and  the  original 
order  of  the  data  is  therefore  lost. 

2.4b  Output.  After  the  user  chooses  to  output  data 
from  IGPS  onto  a  peripheral  device,  he  is  prompted  for  the  type 
of  data  he  wants  to  save,  either  y-only  or  x  &  y.  In  the  first 
case,  only  the  y-axis  values  will  be  output,  whereas  in  the 
second,  both  x-axis  and  y-axis  data  are  transmitted. 

IGPS  then  prompts  for  the  device  the  user  wishes  to 
output  the  data  to.  Possible  responses  are:  one,  the  4907;  two 
the  4924;  and  three,  the  4051  Internal  Tape.  If  the  user  chooses 
the  4907,  IGPS  will  then  ask  for  the  file  names  on  which  to  write 
the  data;  if  he  chooses  the  4924  or  the  internal  tape,  it  asks 
for  the  file  numbers. 

When  writing  to  the  4907,  IGPS  will  create  files  of  the 
appropriate  size  for  the  data  to  be  written;  for  tape  devices, 
however,  the  user  must  use  the  MARK  command  to  create  the  files 
on  his  tape  before  saving  data.  The  tape  files  must  be  s  bytes 
long,  where 

!1 9  +  19  •  m  ,  for  y-only  data 

19  +  37  •  m  ,  for  x  &  y  data 

and  m  is  the  number  of  points  in  the  curve  to  be  written  to  the 
file.  In  practice,  the  default  minimum  file  size  of  768  bytes  is 
adequate  for  most  curves. 

If  the  user  attempts  to  output  data  without  ever  having 
input  data  during  the  session,  an  error  message  is  printed. 
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2.4c  File  Structure.  Data  files  to  be  read  by  IGPS  or 
written  by  IGPS  are  in  ASCII  format.  The  structure  of  an  IGPS- 
compatible  file  is  as  follows:  (1)  the1  first  record  in  the  file 

contains  the  number  of  points  in  the  file;  (2)  each  subsequent 
record  consists  of  either  a  single  number  or  of  a  data  pair, 
depending  on  whether  the  file  is  a  y-only  or  an  x  &  y  data  file. 

2.4d  Address  Number  of  the  4924.  IGPS  assumes  that  the 
4924  Auxiliary  Tape  Drive  is  at  address  2  of  the  General  Purpose 
Interface  Bus  (GPIB).  If  the  user's  4924  is  at  a  different 
address,  he  will  need  to  change  IGPS.  An  example  will  illustrate 
the  required  change:  if  the  address  wanted  is,  say,  4  then  the 

user  should  type  in 


9350  DO  =  4 


This  must  be  done  before  initialization  of  IGPS. 

2.5  User  Key  Number  5:  Curve  Fits. 

This  key  is  used  to  do  least-squares,  linear  fits  to  the 
user's  data.  There  are  several  types  of  fit  available:  (1) 

unconstrained,  where  IGPS  calculates  both  the  slope  and  intercept 
of  the  line  which  best  fits  the  user's  data;  (2)  fixed-slope, 
where  IGPS  calculates  the  intercept  of  the  best  line  with  user- 
defined  slope;  and  (3)  f ixed- intercept ,  where  IGPS  calculates  the 
slope  of  the  best  line  with  user-defined  intercept.  A  fourth 
option  allows  the  user  to  check  the  goodness  of  fit  of  any  line 
he  wishes.  In  all  cases,  IGPS  operates  not  on  the  data  the  user 
inputs  into  IGPS  but  on  the  data  after  it  has  been  transformed 
via  the  two  functions  which  define  the  type  of  paper  the  user 
has  asked  for  (see  section  '.3).  For  instance,  if  the  user 
asked  for  xlog  paper,  the  curve  fit  is  to  the  logarithm  of  the 
x-axis  data  values. 


The  following  equations  are  used  by  IGPS  for  the  different 
curve  fits.  Assume  we  are  given  a  set  of  transformed  data  (xj,y^), 
i  =  l,...,n,  where  n  is  the  number  of  points  in  the  set,  to  which 
we  wish  to  fit  a  line  of  the  form  y  =  ax  +  b.  The  sum  of  the 
squares  of  the  vertical  deviations  between  the  line  and  the  data 
points  is 


D 


Z 


lYi 


(  ax^  +  b)  ]  2  . 


In  order  to  do  an  unconstrained  curve  fit,  in  other  words,  to 
calculate  both  a  and  b,  we  differentiate  D  with  respect  to  a  and 
b  and  set  the  results  to  zero: 
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In  order  to  do  a  fixed-slope  curve  fit,  in  other  words,  in  order 
to  calculate  b  given  a,  we  differentiate  D  with  respect  to  b  and 
set  the  result  to  zero,  which  yields  the  equation  for  b  given 
above.  In  order  to  do  a  f ixed- intercept  fit,  in  other  words,  in 
order  to  calculate  a  given  b,  we  differentiate  D  with  respect  to 
a  and  set  the  result  to  zero,  which  produces: 


a 


xiYi 


b  E 


n 


i=l 


In  all  cases,  IGPS  calculates  and 
goodness  of  fit  of  the  particular 
The  measure  is  known  as  R2  and  is 


displays  a  measure  of 
line  through  the  data 
defined  as: 


the 

points . 
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where  D  is  as  before  and  y  is  the  average  of  y^  ,  i=l,...,n: 

-  1  £ 

y  =  —  Z— *  yi  . 

n  i=l 

From  the  definition  of  ,  it  is  apparent  that  the  closer  r2  is  to 
100,  the  better  the  fit. 

The  prompts  associated  with  user  key  number  5  are 
straightforward  and  require  no  discussion. 

2.6  User  Key  Number  6;  4051. 

This  key  is  used  to  plot  the  current  graph  on  the  screen 
of  the  Tektronix  4051.  If  the  graph  is  not  complete,  in  other 
words,  if  some  item  of  information  needed  to  construct  the  graph 
has  not  been  supplied  to  IGPS,  the  program  will  prompt  the  user 
for  the  missing  item. 

This  key  is  particularly  useful  in  previewing  a  plot: 
the  user  can  look  at  the  plot  on  the  screen  by  depressing  this 
key;  if  he  doesn't  like  what  he  sees,  he  can  use  the  other  keys 
(or  the  shortcuts  described  in  Section  3)  to  alter  the  parameters 
of  the  plot;  finally,  he  can  depress  key  number  6  again  to  quickly 
see  the  results  of  any  changes  he  may  have  made.  Of  course,  if 
the  user  takes  no  action  between  depressions  of  the  key,  then  the 
plot  will  be  re-drawn  on  the  screen  without  changes. 

Once  the  user  is  satisfied  with  the  plot,  he  can  get  a 
copy  of  it  on  the  Tektronix  4662  x-y  plotter  by  depressing  user 
key  number  7 . 

2.7  User  Key  Number  7:  4662. 

This  key  is  used  to  obtain  a  hard  copy  of  the  current 
graph  on  Tektronix  4662  x-y  plotter.  Before  using  this  key,  the 
user  should  prepare  the  4662  for  plotting  by  placing  a  piece  of 
paper  in  the  lower  left-hand  corner  of  the  plotting  area  of  the 
4662.  If  the  piece  of  paper  is  11  inches  long  and  eight  and  one- 
half  inches  high,  and  if  the  hardware  limits  of  the  4662  are  left 
at  their  default  values,  the  plot  produced  by  IGPS  will  be  centered 
on  the  piece  of  paper  (see  Figure  1.1).  To  produce  plots  of  dif¬ 
ferent  sizes,  the  user  can  either  change  the  hardware  limits  of 
the  4662  or  he  can  modifiy  the  definition  of  the  viewport  in  IGPS 
(see  the  listing  of  IGPS  in  Appendix  A  at  lines  6620  through  6660). 
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IGPS  assumes  that  the  4662  plotter  is  at  address  1  of 
the  GPIB.  If  the  user's  4662  is  at  a  different  address,  he  will 
need  to  change  line  28  of  the  program  to  read 

-  28  D= 

where  the  new  address  goes  to  the  right  of  the  equal  sign. 

2.8  User  Key  Number  8;  Storage  and  Retrieval. 

This  key  is  used  to  store  onto  a  peripheral  or  retrieve 
from  a  peripheral  all  the  information  required  to  construct  a 
particular  plot.  This  feature  is  useful  in  that  it  allows  a 
complex  plot  to  be  saved  during  one  session  with  IGPS  and  then 
quickly  reconstructed  at  a  later  session  without  the  user  having 
to  re-input  all  the  parameters  which  make  up  the  plot;  instead, 
the  entire  plot  can  be  recreated  by  just  depressing  this  one  key 
and  responding  to  a  few  prompts. 

After  the  user  depresses  key  number  8,  IGPS  will  prompt 
him  for  the  type  of  operation  to  carry  out  (either  storage  or 
retrieval).  The  user  is  then  prompted  for  the  device  he  wishes 
to  use.  There  are  three  possible  responses:  one,  meaning  the 
4907  file  manager;  two,  meaning  the  4924  Tape  Drive;  and  three, 
meaning  the  internal  drive.  If  the  user  selected  the  4907,  he  is 
then  prompted  for  the  name  of  the  file  to  use;  if  he  selected  the 
4924  or  the  interval  tape  drive,  he  is  prompted  for  the  number  of 
the  file. 


If  the  user  attempts  to  store  an  incomplete  plot,  an  error 
message  is  printed  and  no  further  action  is  taken. 

When  storing  a  plot  on  the  4907,  IGPS  will  create  a 
binary  file  of  the  appropriate  size  for  the  data  to  be  written; 
for  tape  devices,  the  user  must  use  the  MARK  command  to  create 
the  file  on  his  tape.  The  tape  file  must  be  S  bytes  long,  where 
S  is  approximately  given  by  the  expression 

S  =  1000  +  16  *  n  *  mQ, 

where  n  is  the  number  of  data  curves  in  the  plot  and  m0  is  the 
number  of  points  in  the  curve  having  the  greatest  number  of 
points;  for  instance,  storing  Figure  1.1  requires  approximately 
2000  bytes. 

The  remarks  in  2.4d  about  the  address  of  the  4924  apply 
here  as  well. 

2.9  User  Key  Number  9:  Data  Generation. 

This  key  can  be  used  to  generate  x  axis  and  y  axis  data 
from  user-defined  parametric  functions  rather  than  inputting  them 
from  a  peripheral  or  the  keyboard. 
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After  the  user  depresses  this  key,  IGPS  will  prompt 

with: 

DATA  TYPE  (1=Y  ONLY,  2=X&Y): 

NUMBER  OF  CURVES: 

NUMBER  OF  POINTS  PER  CURVE  (AT  LEAST  2): 

MIN  AND  MAX  OF  PARAMETRIC  VARIABLE: 

The  meaning  of  the  DATA  TYPE  prompt  is  as  before:  a 
response  of  1  is  a  request  to  IGPS  to  sort  the  y  axis  data  in 
increasing  order  and  to  generate  x  axis  data  as  cumulative 
fractions. 

IGPS  will,  after  the  user  responds  to  the  prompts,  create 
the  data  as 


xi,j 

II 

Tl 

X 

t  tG  +  s  •  (  j 

-  1)] 
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Yi,j 
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tQ  and  tj  are  the  minimum  and  maximum  of  the  parametric  variable, 
and  Fx  and  Fy  are  the  user-definable  functions.  These  two 
functions  must  be  defined  by  the  user  before  he  depresses  user 
key  number  9;  he  must  do  so  by  replacing  lines  3420  and  3430  in 
the  program  with  his  own  functions:  line  3240  is  for  Fx  and  must 
be  of  the  form 

3420  DEF  FNX ( T ) = 

while  line  3430  is  for  Fy  and  must  be  of  the  form 

3430  DEF  FNY ( T ) = 

By  default,  the  definitions  are 

3420  DEF  FNX ( T )  =  RO ( I ) *  SIN ( T ) 

3430  DEF  FNY ( T ) =RO (I ) *SIN ( T ) *COS ( T ) 
where  RO  is  an  array  with  five  entries  such  that 

ROi  =  i ,  i=l , . . . , 5 . 
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Therefore,  if  the  user  does  not  alter  the  definitions,  and  if  he 
uses  the  following  sequence 

DATA  TYPE  (1=Y  ONLY,  2=X&Y):  2 

NUMBER  OF  CURVES:  5 

NUMBER  OF  POINTS  PER  CURVE:  21 

MIN  AND  MAX  OF  INDEPENDENT  VARIABLE:  0,  6.28 

then,  with  appropriate  choices  of  limits,  labels,  curve  types  and 
so  on,  he  should  be  able  to  reproduce  Figure  2.15. 


Figure  2.15  Sample  Plot  Using  Parametric  Functions 
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2.10  User  Key  Number  10:  Menu. 


Whenever  this  key 
on  the  screen  of  the  4051; 


is  depressed,  IGPS  will  plot  the  menu 
see  Figure  2.1. 


Next  page  is  blank 


41 


3 .  SHORTCUTS 


Suppose  you  are  previewing  a  plot  on  the  screen  of  the 
Tektronix  4051  and  you  decide  to  change  tne  type  of  curve  marker 
being  used  with  one  of  the  curves.  One  way  to  do  this  is  to 
depress  user  key  number  4  and  to  go  through  the  input  procedure; 
this,  however,  is  clumsy  and  can  be  avoided  since  the  variable 
controlling  the  type  of  curve  marker  is  directly  available  to  the 
user:  it  is  in  fact  the  array  LO,  with  as  many  entries  as  there 

are  curves,  and  with  the  i^  entry  of  LO  being  an  integer  between 
one  and  nine  (see  Figure  2.14)  corresponding  to  the  type  of  curve 
marker  to  be  used  with  the  i^“  curve.  Therefore,  to  change  the 
type  of  curve  marker  for  a  given  curve,  the  user  need  only  change 
the  entry  in  LO  corresponding  to  that  curve.  For  instance,  to 
draw  the  third  curve  as  a  solid  line,  the  user  should  type  in 
LO  (  3  )  =  4 . 


The  same  procedure  can  be  used  with  any  of  the  features 
of  the  plot  being  previewed.  To  find  out  the  name  of  the  variable 
which  controls  the  feature  of  the  plot  which  you  want  to  change, 
consult  Table  3-4. 


TABLE  3-1  DICTIONARY  OF  VARIABLES 


Name 

Type 

Description 

Nl 

S imple 

Number  of 

entries  in  XI 

X1(N1) 

Array 

Locations 

of  x  axis  tick 

marks 

N3 

Simple 

Number  of 

entries  in  Y1 

Y1(N3) 

Array 

Locations 

of  y  axis  tick 

marks 

N2 

S imple 

Number  of 

entries  in  A$ 

A$a 

String 

X  axis  tick  mark  labels 

N4 

Simple 

Number  of 

entries  in  B$ 

B$a 

String 

Y  axis  tick  mark  labels 

x$ 

String 

X  axis  label 

Y$ 

String 

Y  axis  label 

z$ 

String 

Title  of 

plot 

a  see  remarks 
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TABLE  3-1  DICTIONARY  OF  VARIABLES  (CONTINUED) 


Name 

TYPe 

Description 

AO 

Simple 

Frame  type: 

0  -  None 

1  -  Axes 

2  -  Grid 

3  -  Border 

FO 

Simple 

Color  for  frame  and  labels 
(must  be  integer  between 
one  and  eight) 

SO 

Simple 

Pen  speed  (must  be  number 
between  10  and  570) 

Pa 

Simple 

Paper  type: 

1  -  linear 

2  -  Xlog 

3  -  Ylog 

4  -  Loglog 

5  -  User-defined 

6  -  Uniform  probability 

7  -  Normal  probability 

8  -  Lognormal  probability 

9  -  Rayleigh  probability 
10  -  Weibull  probability 

N 

Simple 

Number  of  curves 

MO 

Simple 

Number  of  points  in  the 
curve  with  the  greatest 
number 

M  (N ) 

Array 

M(I)  contains  the  number  of 
points  in  the  i^h  curve 

X (N ,MO ) a 

Array 

X(I,J)  contains  the  jth 
data  point  in  the  i^n 
curve 

Y ( N ,MO ) a 

Array 

Y(I,J)  contains  the  jth 
data  point  in  the  itn 

curve 


see 


remarks 
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TABLE  3-1  DICTIONARY  OF  VARIABLES  (CONTINUED) 


Name 

TYPe 

Description 

LO  ( N ) 

Array 

LO(I)  contains  the  curve 
type  for  the  i^h  curve: 

1  -  Dashed  line 

2  -  Dashed-dotted  line 

3  -  Dotted  line 

4  -  Solid  line 

5  -  Crosses 

6  -  X-shaped  markers 

7  -  Squares 

8  -  Triangles 

9  -  Upside-down  triangles 

F$a 

String 

Legends  for  the  N  curves 

CO  ( N ) 

Array 

CO(I)  contains  the  color  of 
the  i^n  curve  and  legend  (musl 
be  integer  between  one  and 
eight) 

C 

Simple 

Curve  fit: 

0  -  No 

7^0  -  Yes 

A  (N  ) a 

Array 

A(  I )  contains  the  slope  of 
the  line  fit  to  the  i^h 
set  of  data 

B(N)a 

Array 

B ( I )  contains  the  intercept 
of  the  line  fit  to  the  ita 
set  of  data 

R(N)a 

-a - ; — 

Array 

R(I)  contains  the  r2  - 
value  associated  with  the 
fit  to  the  i^h  set  of 
data 

a  see  remarks 


REMARKS : 

1.  The  strings  A$  and  B$  contain  the  tick  mark  labels  in  compact 
form,  which  is  to  say,  A$  consists  of  N2  labels  juxtaposed  and 
separated  by  pound  signs  (#).  Therefore,  the  tick  mark  labels 
associated  with  Figure  2.2  could  have  been  created  by  typing  in 
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N2  =  5 


A$="0.00#0.25#0.05#0.75#1.00#" 

2.  Changes  to  the  variable  P  must  be  followed  by  the  statement 
RUN  4410. 

3.  Unused  entries  of  the  arrays  X  and  Y  are  zero-filled. 

4.  F$  contains  the  N  legends  in  compact  form  (see  remark  1  for  a 
definition  of  compact).  The  legends  in  Figure  1.1  could  have 
been  created  by  the  statement 

F$= " MARK  I#  MARK  II#" 

5.  The  contents  of  the  arrays  A,  B,  and  R  are  only  meaningful 
if  a  curve-fit  has  been  performed.  Otherwise,  they  are  zero- 
filled. 
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4. 


POSTSCRIPT 


The  author  has  attempted  to  make  IGPS  as  useful  and 
versatile  as  possible  and  he  welcomes  recommendations  for 
improving  it. 


Next  page  is  blank 
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APPENDIX  A 
PROGRAM  LISTING 


SET  UP  USER-DEFINABLE  KEYS  AND  PROGRAM  INITIALIZATION 
(NOTE:  4662  ASSUMED  AT  ADDRESS  1  OF  GPIB;  SEE  LINE  28) 


1 

GOSUB 

9330 

2 

GO 

TO 

3640 

4 

GO 

TO 

100 

8 

GO 

TO 

1310 

12 

GO 

TO 

4150 

16 

GO 

TO 

1410 

20 

GO 

TO 

5560 

24 

D=32 

25 

GO 

TO 

6390 

28 

D=1 

29 

GO 

TO 

6390 

32 

GO 

TO 

4920 

36 

GO 

TO 

3370 

40 

GO 

TO 

3640 

USER  KEY  Is  SET  UP  AXIS  LIMITS  AND  LABELS 


100  PAGE 

110  DELETE  XI , Yl 
120  A$  = " " 

130  PRINT  "UNIFORM  TICKS  ON  THE  X-AXIS  ( 1= YES , 2=NO) : " ; 
140  INPUT  L 

SELECT  TYPE  OF  X-AXIS  LABELING 
150  GO  TO  L  OF  160,400 

1. AUTOMATIC  X  AXIS 

160  REM 
170  Ml=0 

180  PRINT  "XMIN : " ; 

190  GOSUB  1010 

200  D1=VAL(W$) 

210  PRINT  "XMAX:"; 

220  GOSUB  1010 

230  D2=VAL( W$ ) 

240  PRINT  "MAJOR  TICKS:"; 

250  GOSUB  1010 

260  D4=VAL(W$) 

270  PRINT  "MINOR  TICKS:"; 

280  INPUT  D3 

290  Nl=l+INT( ( D2-D1 )/D3 ) 

300  DIM  X1(N1) 

310  FOR  1=1  TO  Nl 

320  X1(I)=D1+(I-1)*D3 
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330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 


NEXT  I 

N2=1+INT( (D2-D1 )/D4 ) 

FOR  I =Dl  TO  D2  STEP  D4 
GOSUB  1060 
A$=A$  &W$ 

NEXT  I 
GO  TO  540 

2.  USER-CONTROLLED  X  AXIS 


REM 

PRINT  "NUMBER  OF  TICK  MARKS  (AT  LEAST  2):"; 

INPUT  Nl 
DIM  XI (Nl ) 

PRINT  "LOCATION  OF  TICK  MARKS:"; 

INPUT  XI 

PRINT  "NUMBER  OF  TICK  MARK  LABELS:"; 

INPUT  N2 

IF  N2=0  THEN  530 

PRINT  "LOCATION  OF  TICK  MARK  LABELS:"; 

INPUT  W$ 

GOSUB  1230 
A$=W  $ 

REM 
REM 
B$=  "  " 

PRINT  "UNIFORM  TICKS  ON  THE  Y-AXIS  ( 1 =YES , 2=NO ) : " ; 
INPUT  L 


SELECT  TYPE  OF  Y-AXIS  LABELING 


GO  TO  L  OF  590,830 


REM 


AUTOMATIC  Y  AXIS 


M1  =  0 

PRINT  " YMIN : " ; 

GOSUB  1010 
Dl =VAL ( W$ ) 

PRINT  "YMAX: " ; 

GOSUB  1010 
D2=VAL ( W$ ) 

PRINT  "MAJOR  TICKS:"; 
GOSUB  1010 
D4  =VAL ( W $ ) 

PRINT  "MINOR  TICKS :" ; 
INPUT  D3 

N3=l +INT ( ( D2-D1 ) /D3 ) 
DIM  Y1 ( N3 ) 

FOR  1=1  TO  N3 

Y1 ( I ) =D1  +  ( I  — 1 ) *D  3 
NEXT  I 

N4=1+INT ( (D2-D1 )/D4 ) 
FOR  I=D1  TO  D2  STEP  D4 
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790 

GOSUB 

1060 

800 

B$  =  B$  &W$ 

810 

NEXT  I 

820 

GO  TO  970 

2. 

USER-CONTROLLED  Y 

AXIS 

830 

REM 

840 

PRINT  "NUMBER 

OF  TICK  MARKS  (AT 

LEAST 

850 

INPUT  N3 

860 

DIM  Yl(N3) 

870 

PRINT  "LOCATION  OF  TICK  MARKS:" 

/ 

880 

INPUT  Y1 

890 

PRINT  "NUMBER 

OF  TICK  MARK  LABELS:"; 

900 

INPUT  N4 

910 

IF  N4=0 

THEN 

960 

920 

PRINT 

"LOCATION  OF  TICK  MARK 

LABELS 

930 

INPUT 

W$ 

940 

GOSUB 

1230 

950 

B$  =  W  $ 

960 

REM 

970 

REM 

980 

K  (  1 )  =  1 

990 

PAGE 

1000 

RETURN 

ROUTINE  TO  CALCULATE  THE  NUMBER  OF  DIGITS  TO  THE  RIGHT  OF 
DECIMAL  POINT 

1010  INPUT  W $ 

1020  IF  POS( W$ , " . " , 1 )=0  THEN  1040 
1030  Ml=Ml  MAX  LEN(W$)-POS(W$,".",l) 

1040  REM 
1050  RETURN 

ROUTINE  TO  CONSTRUCT  TICK  MARK  LABEL  WITH  A  GIVEN  NUMBER 
PLACES,  Ml,  TO  THE  RIGHT  OF  THE  DECIMAL  POINT 

1060  W$=STR( ABS ( I ) * 10“M1 ) 

1070  W$=SEG(W$,2,LEN(W$)-1) 

1080  IF  Ml=0  THEN  1170 

1090  IF  Ml <=  LEN ( W$ )  THEN  1120 

1100  U$=SEG( "00000000000" ,1 ,Ml-LEN( W$ ) ) 

1110  W$=U$&W$ 

1120  REM 

1130  W$=REP ( " . " , LEN ( W$ ) -Ml+1 , 0 ) 

1140  IF  ABS ( I ) => 1  THEN  1160 
1150  W$=REP( "0" ,1 ,0) 

1160  REM 
1170  REM 

1180  IF  I=>0  THEN  1200 
1190  W$=REP( ,1 ,0 ) 

1200  REM 
1210  W$=W$S,"#" 
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1220  RETURN 


ROUTINE  TO  REPLACE  BLANKS  OR  COMMAS  WITH  POUND  SIGNS  (#) 

1230  FOR  1=1  TO  LEN ( W$ ) 

1240  U$=SEG(W$, 1,1 ) 

1250  IF  U$<>"  "  AND  U$<>"  ,"  THEN  1270 
1260  W$=REP("#" ,1,1) 

1270  REM 
1280  NEXT  I 
1290  W$=W$&"#" 

1300  RETURN 

USER  KEY  2:  INPUT  TITLES 


1310  PAGE 

1320  PRINT  "X-AXIS  LABEL:"; 

1330  INPUT  X$ 

1340  PRINT  "Y-AXIS  LABEL:"; 

1350  INPUT  Y$ 

1360  PRINT  "PLOT  TITLE:"; 

1370  INPUT  Z$ 

1380  K ( 2 ) =1 
1390  PAGE 
1400  RETURN 

USER  KEY  4:  INPUT  AND  OUTPUT  OF  DATA 


1410  PAGE 

1420  PRINT  "INPUT  (1)  OR  SAVE  (2):"; 

1430  INPUT  L 

1440  PRINT  "DATA  TYPE  ( 1=Y  ONLY , 2=X&Y ) : " ; 

1450  INPUT  NO 

BRANCH  ON  WHETHER  INPUT  OR  OUTPUT 

1460  GO  TO  L  OF  1470,2110 
1470  REM 

1480  PRINT  "UNIT  NUMBER  ( 1=4907 , 2=4924 , 3=INTERNAL  TAPE , 4=KEYBOARD ) 
1490  INPUT  L 

1500  PRINT  "NUMBER  OF  CURVES:"; 

1510  INPUT  N 

1520  PRINT  "MAXIMUM  NUMBER  OF  POINTS  PER  CURVE:"; 

1530  INPUT  M0 

1540  DELETE  X,Y,M 

1550  DIM  X ( N , M0 ) , Y ( N , M0 ) , M ( N ) 

1560  X=0 

1570  Y=0 

1580  FOR  1=1  TO  N 
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1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1  880 

1890 

1900 

1910 

1920 

1930 

1  940 

1950 

1960 

1970 

1980 

1  990 

2000 


BRANCH  ON  DEVICE  TYPE 


GO  TO  1  * ( L=1 ) +2 * ( L=2  OR  L=3)+3*(L=4)  OF  16  00,1760,1920 

1.  INPUT  FROM  4907  FILE  MANAGER 


REM 

PRINT  "FILE  NAME  FOR  CURVE  NUMBER  " ; I ; " : " ; 

INPUT  W$ 

OPEN  W$ ; 1 , "R" ,U$ 

INPUT  #1 : M ( I ) 

FOR  J=1  TO  M ( I ) 

GO  TO  NO  OF  1670,1700 
REM 

INPUT  #1 : Y( I , J ) 

GO  TO  1720 
REM 

INPUT  #1 :X( I , J ) , Y ( I , J ) 

REM 

NEXT  J 
CLOSE  1 
GO  TO  2070 

2.  INPUT  FROM  4924  TAPE  DRIVE  OR  INTERNAL  DRIVE 
(NOTE:  4924  ASSUMED  AT  ADDRESS  2  OF  GPIB;  SEE 

LINE  9350) 

REM 

PRINT  "FILE  NUMBER  FOR  CURVE  NUMBER  " ; I ;  "  :  "  ; 

INPUT  Dl 

D=D0* ( L=2 )  +  33  * ( L=3 ) 

FIND  @D: Dl 
INPUT  @D:M( I ) 

FOR  J=1  TO  M ( I ) 

GO  TO  NO  OF  1840,1870 
REM 

INPUT  @D: Y ( I , J ) 

GO  TO  1890 
REM 

INPUT  @D: X ( I , J ) , Y ( I , J ) 

REM 

NEXT  J 
GO  TO  2070 


3.  INPUT  FROM  KEYBOARD 


REM 

PAGE 

PRINT  "NUMBER  OF  POINTS  IN  CURVE  NUMBER  " ; I ;  "  :  " ; 
INPUT  M ( I ) 

FOR  J  =  1  TO  M (  I ) 

GO  TO  NO  OF  1980,2020 
REM 

PRINT  " Y ( " ; I ; " , "; J;  ") :  "; 

INPUT  Y ( I , J ) 
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2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 


GO  TO  2050 
REM 

PRINT  "X( ";I; " , " ;J; " ) ,Y( " ; I; " , " ; J; " ) : " ; 

INPUT  X(I,J) ,Y(I ,J) 

REM 

NEXT  J 
REM 

NEXT  I 
GOSUB  2650 
GO  TO  2630 
REM 

IF  K ( 4  )  =  0  THEN  2570 

PRINT  "UNIT  NUMBER  ( 1=4907 , 2=4924 , 3=INTERNAL  TAPE):"; 

INPUT  L 

FOR  1=1  TO  N 

BRANCH  ON  DEVICE  TYPE 

GO  TO  1* ( L=1 )  +  2 * ( L=2  OR  L=3)  OF  2170,2350 

1.  OUTPUT  TO  4907  FILE  MANAGER 


REM 

PRINT  "FILE  NAME  FOR  CURVE  NUMBER  " ; I ; " : " ; 

INPUT  W$ 

KILL  W$ 

CREATE  W$,"A";19+(19*(N0=1)+37*(N0=2) )*M(I) ,0 
OPEN  W$ ; 1 , "F" ,U$ 

PRINT  # 1 :M ( I ) 

FOR  J=1  TO  M ( I ) 

GO  TO  NO  OF  2260,2290 
REM 

PRINT  # 1 : Y ( I , J ) 

GO  TO  2310 
REM 

PRINT  # 1 : X ( I , J ) , Y( I , J ) 

REM 

NEXT  J 
CLOSE  1 
GO  TO  2530 

2.  OUTPUT  TO  4924  TAPE  DRIVE  OR  INTERNAL  DRIVE 

(NOTE:  4924  ASSUMED  AT  ADDRESS  2  OF  GPIB;  SEE 

LINE  9350) 

REM 

PRINT  "FILE  NUMBER  FOR  CURVE  NUMBER  "  ;  I ;  "  :  "  ; 

INPUT  Dl 

D=D0*(L=2)+33*(L=3) 

KILL  @D:D1 
FIND  @D : Dl 
PRINT  @D:M( I ) 

FOR  J=1  TO  M ( I ) 
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2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 


2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 


2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 


2870 

2880 


GO  TO  NO  OF  2440,2470 
REM 

PRINT  @D: Y( I , J ) 

GO  TO  2490 
REM 

PRINT  @D:X( I , J ) , Y ( I , J ) 

REM 

NEXT  J 
PRINT  @D,  2 : 

CLOSE 
REM 
NEXT  I 
PAGE 

GO  TO  2620 

ERROR  MESSAGE  IF  OUTPUT  ATTEMPTED  WITHOUT  INPUT 


REM 

PAGE 

PRINT  @32,21:35,50 

PRINT  "NO  INPUT  YET — NOTHING  TO  SAVE" 
HOME 
REM 
REM 

RETURN 

INPUT  CURVE  TYPES 

PAGE 

PRINT  "CURVE  TYPES  AVAILABLE:" 

PRINT 

PRINT  "  1. DASHED  LINE" 

PRINT  "  2. DASHED- DOTTED  LINE" 

PRINT  "  3. DOTTED  LINE" 

PRINT  "  4. SOLID  LINE" 

PRINT  "  5. CROSSES  " 

PRINT  "  6 . X 1 S" 

PRINT  "  7. SQUARES  " 

PRINT  "  8. TRIANGLES  " 

PRINT  "  9. FLASHES" 

PRINT 

DELETE  L0 , CO , A, B , R 

DIM  L0(N),C0(N) ,A(N) ,B(N) ,R(N) 

A=0 

B=0 

C=0 

FOR  1=1  TO  N 

PRINT  "TYPE  FOR  CURVE  NUMBER  "  ;  I ;  "  :  " ; 
INPUT  L0 ( I ) 

NEXT  I 

INPUT  LEGENDS 

PAGE 
F$="  " 


59 


2890 

2900 

2910 

2920 

2930 

2940 


2950 

2960 

2970 

2980 

2990 


3000 

3010 

3020 

3030 

3040 

3050 


3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 


FOR  1=1  TO  N 

PRINT  "LEGEND  FOR  CURVE  NUMBER  "  ;  I  ;  "  :  "  ; 
INPUT  W$ 

F$=F$&W$ 

F$=F$& "# " 

NEXT  I 

INPUT  COLORS 

PAGE 

FOR  1=1  TO  N 

PRINT  "COLOR  FOR  CURVE  NUMBER  "  ;  I ;  "  :  "  ; 
INPUT  CO (I) 

NEXT  I 

GENERATE  X-AXIS  DATA  IF  REQUESTED 

IF  N0=  2  THEN  3330 
FOR  1=1  TO  N 

FOR  J=1  TO  M( I ) 

X(I, J)=100*J/(1+M(I) ) 

NEXT  J 
NEXT  I 


SORT  Y-AXIS  DATA  IF  REQUESTED  (USES  BUBBLE  SORT) 

FOR  1=1  TO  N 

IF  M ( I ) <= 1  THEN  3310 
FOR  J=  2  TO  M ( I ) 

IF  Y( I ,J-1 ) <=Y(I ,J)  THEN  3290 
Dl=Y (  I ,  J  ) 

L=  J  +  l 
D2  =  l 
REM 

L=L- 1 

IF  L=1  THEN  3240 

IF  Y  (  I , L— 1 ) <=Dl  THEN  3190 
Y(I,L)=Y(I,L-1) 

GO  TO  3220 
REM 

Y (  I  , L ) =Dl 
D2=0 
REM 

GO  TO  3270 
REM 

Y(I,1)=D1 

D2=0 

REM 

IF  D2  THEN  3130 
REM 

NEXT  J 
REM 

NEXT  I 
REM 
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3340  K ( 4 ) =1 
3350  PAGE 
3360  RETURN 

USER  KEY  9:  GENERATE  DATA  FROM  FUNCTIONS 


3370  PAGE 
3380  DIM  R0 ( 5 ) 

3390  RESTORE  3400 
3400  DATA  1,2, 3,4, 5 
3410  READ  R0 

DEFINE  FUNCTIONS 

3420  DEF  FNX(T)=R0(I)*SIN(T) 

3430  DEF  FNY(T)=R0(l)*SIN(T) *COS ( T ) 

3440  PRINT  "DATA  TYPE  ( 1=Y  ONLY , 2=X&Y ) : " ; 

3450  INPUT  NO 

3460  PRINT  "NUMBER  OF  CURVES:"; 

3470  INPUT  N 

3480  PRINT  "NUMBER  OF  POINTS  PER  CURVE  (AT  LEAST  2):"; 

3490  INPUT  M0 

3500  DELETE  X,Y,M 

3510  DIM  X ( N , M0 ) , Y ( N , M0 ) ,M(N) 

3520  M=M0 

3530  PRINT  "MIN  AND  MAX  OF  PARAMETRIC  VARIABLE:"; 

3540  INPUT  TO , Tl 
3550  S=(T1-T0)/(M0-1 ) 

GENERATE  DATA  (NOTE:  I  INDEXES  THE  CURVES,  J  THE  POINTS) 

3560  FOR  1=1  TO  N 

3570  FOR  J=1  TO  M0 

3580  X( I, J)=FNX(T0+S* (J-l ) ) 

3590  Y(I ,J)=FNY(T0+S*( J-l) ) 

3600  NEXT  J 
3610  NEXT  I 
3620  GOSUB  2650 
3630  RETURN 

USER  KEY  10:  DRAW  MENU 


3640  D=3 2 
3650  PAGE  @D : 

3660  DIM  K$ (80) 

3670  K$= "  LIMITS  LABELS  FORMAT  IO  FITS  " 

3680  K$  =  K$ & "  4051  4662  S&R  FUNCTION  MENU 

3690  L=6 0 

3700  D2=L 

3710  FOR  J=1  TO  2 

3720  Dl=-4 

3730  FOR  1=1  TO  5 
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3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 


Dl=Dl+20 

PRINT  @D,21:D1,D2 

PRINT  @D , 20 : Dl+1 6 , D2 ,D1+16 , D2-4 , Dl , D2-4 , Dl , D2 
NEXT  I 
D2=D2-17 . 6 


NEXT  J 
PRINT 

@D,21:23.8,L-8.4 

PRINT 

@D : "1" 

PRINT 

@D,21:43.8,L-8.4 

PRINT 

@D: "2" 

PRINT 

@D,21:63.8,L-8.4 

PRINT 

@D: "3" 

PRINT 

@D,21:83.8,L-8.4 

PRINT 

@D: "4" 

PRINT 

@D, 21:103. 8, L- 8. 4 

PRINT 

@D: "5" 

PRINT 

@D,21:16,L-10.8 

PRINT 

@D,20: 50 ,L-10.8 

PRINT 

@D,  21 : 52 ,L-12 

PRINT 

@D : "USER  DEFINABLE 

PRINT 

@D,21:78,L-10.8 

PRINT 

@D,20: 112, L- 10. 8 

PRINT 

@D,21:23.8,L-26 

PRINT 

@D: "6" 

PRINT 

@D, 21:43. 8, L-26 

PRINT 

@D: "7" 

PRINT 

@D, 21: 63.8, L-26 

PRINT 

@D: "8" 

PRINT 

@D, 21:83. 8, L-26 

PRINT 

@D: "9" 

PRINT 

@D, 21:102. 5, L-26 

PRINT 

@D : "10" 

FOR  J= 

=  1  TO  2 

FOR 

1=1  TO  5 

PRINT  @D, 21: 17+20* (1-1) , L-3 . 2-17 . 6* ( J-l ) 
U$=SEG(K$, (J-l )* 40+ (1-1 ) *8+1,8) 

PRINT  @D : U$ 

NEXT  I 
NEXT  J 
HOME  @D: 

RETURN 


USER  KEY  3:  FORMAT  DEFINITION 


PAGE 

PRINT  "FRAME  TYPE  ( 0=NONE , 1=AXES , 2=GRID , 3= BORDER) : 
INPUT  AO 

PRINT  "COLOR  FOR  FRAME  AND  LABELS  (1-8):"; 

INPUT  F0 

PRINT  "PEN  SPEED  (10-570  MM/S ) : " ; 

INPUT  SO 
PRINT 

PRINT  "  PAPER  TYPES:  X-AXIS  Y 


«  . 
r 


-AXIS" 
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4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 

4730 

4740 


PRINT 

PRINT 

"  1. 

LINEAR: 

X 

Y" 

PRINT 

"  2. 

XLOG: 

LOG  ( X ) 

Y" 

PRINT 

"  3. 

YLOG : 

X 

LOG 

(Y)" 

PRINT 

"  4. 

LOGLOG: 

LOG ( X ) 

LOG 

(Y)" 

PRINT 

"  5. 

USER: 

•p 

?" 

PRINT 

"  6. 

UNIFORM: 

X/100 

Y" 

PRINT 

"  7. 

NORMAL: 

N ( X/l 00 ) 

Y" 

PRINT 

"  8. 

LOGNORMAL: 

N ( X/100 ) 

LOG 

(Y)" 

PRINT 

"  9. 

RAYLEIGH:  SQR ( -2 *LOG ( l-X/100 ) ) 

Y" 

PRINT 

"10. 

WEIBULL: 

LOG (-LOG (l-X/100 ) ) 

LOG 

(Y)" 

PRINT 

PRINT 

"PAPER 

TYPE: " ; 

INPUT 

P 

GOSUB 

4410 

K ( 3 ) =1 

RETURN 

BRANCH  ON  PAPER 

TYPE  SELECTED  FOR  FUNCTION 

DEFINITION 

GO  TO 

P  OF  4420,4460,4500 

,4540,4580,4620,4660 

,4740,4820,4860 

REM 

DEF 

FNA(X) 

=x 

DEF 

FNB ( Y ) 

=  Y 

GO  TO  4890 

REM 

DEF  FNA ( X )  =  LOG ( X ) 

DEF  FNB ( Y ) =Y 
GO  TO  4890 
REM 

DEF  FNA(X)=X 
DEF  FNB ( Y ) =LOG ( Y ) 

GO  TO  4890 
REM 

DEF  FNA( X ) =LOG ( X ) 

DEF  FNB ( Y ) =LOG ( Y ) 

GO  TO  4890 
REM 

DEF  FNA( X ) =X 
DEF  FNB ( Y ) =Y 
GO  TO  4890 
REM 

DEF  FNA( X ) =X/1 00 
DEF  FNB ( Y ) =Y 
GO  TO  4890 
REM 

DEF  FNA( X ) =FNC ( ( l-X/100 ) * ( X>50 ) +X/100 * ( X<= 50 ) ) * ( ( X> 50 ) - ( X<=50 ) ) 
DEF  FNC(X)=FND(X)-FNE( FND(X) )/FNF( FND(X) ) 

DEF  FND ( X ) =SQR ( LOG ( 1/ ( X*X ) ) ) 

DEF  FNE (X)=2.515517+( 0.802853+0 .010328*X) *X 
DEF  FNF(X)=1+(1 . 432788+ ( 0.189269+0 . 00 13 08 *X ) *X ) *X 
DEF  FNB ( Y ) =Y 
GO  TO  4890 
REM 
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4750 

4760 

4770 

4780 

4790 

4800 

4810 

4820 

4830 

4840 

4850 

4860 

4870 

4880 

4890 

4900 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4980 

4990 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5100 


DEF  FNA( X )=FNC ( ( l-X/100 ) * (X>50 )+X/100* (X<=50 ) )*( ( X>50 )- (X<=50 ) ) 
DEF  FNC (X ) =FND ( X ) -FNE ( FND ( X ) )/FNF( FND(X) ) 

DEF  FND (X ) =SQR( LOG ( 1/ (X*X ) ) ) 

DEF  FNE (X) =2. 5 1551 7+ (0.802853+0.0 10328 *X)*X 

DEF  FNF(X)=1+(1 .432788+ (0.189269+0 .001 308 *X)*X)*X 

DEF  FNB ( Y ) =LOG ( Y ) 

GO  TO  4890 
REM 

DEF  FNA(X)=SQR(-2*LOG( l-X/100 ) ) 

DEF  FNB ( Y ) =Y 
GO  TO  4890 
REM 

DEF  FNA(X)=LOG(-LOG( l-X/100 ) ) 

DEF  FNB ( Y ) =LOG ( Y ) 

REM 

PAGE 

RETURN 


USER  KEY  8:  STORE  AND  RETRIEVE  PLOTS 


PAGE 

PRINT  "STORE  (1)  OR  RETRIEVE  (2):"; 

INPUT  L 

PRINT  "UNIT  NUMBER  ( 1=49 07 , 2=49 24 , 3= INTERNAL  TAPE):"; 

INPUT  Dl 

BRANCH  ON  STORAGE  OR  RETRIEVAL 

GO  TO  L  OF  4980,5310 
REM 

IF  SUM(  K  )  05  THEN  5240 

BRANCH  ON  DEVICE  TYPE 

GO  TO  1* (Dl=l )+2* (Dl=2  OR  Dl=3)  OF  5010,5110 
1.  STORAGE  ON  4907  FILE  MANAGER 

REM 

PRINT  "FILE  NAME:"; 

INPUT  W$ 

KILL  W$ 

CREATE  W$ ; 845+8* (N1+N3+N* ( 2*M0+5 ) ) ,0 
OPEN  W$ ; 1 , " F" ,U$ 

WRITE  #1:N1,N2,N3,N4,A$,B$,X$,Y$,Z$,N,M0,F$,A0,F0,S0,P,C 
WRITE  #1:X1,Y1,X,Y,M,L0,C0,A,B 
CLOSE  1 
GO  TO  5210 
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5110 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5190 

5200 

5210 

5220 

5230 

5240 

5250 

5260 

5270 

5280 

5290 

5300 

5310 

5320 

5330 

5340 

5350 

5360 

5370 

5380 

5390 

5400 

5410 

5420 

5430 

5440 

5450 

5460 


2.  STORAGE  ON  4924  TAPE  DRIVE  OR  INTERNAL  DRIVE 
(NOTE:  4924  ASSUMED  AT  ADDRESS  2  OF  GPIB;  SEE 
LINE  9350) 


REM 

PRINT  "FILE  NUMBER:"; 

INPUT  D2 

D=D0* (Dl=2)+33* (Dl=3) 

KILL  @D:D2 
FIND  @D:D2 

WRITE  @D:Nl,N2,N3,N4,A$,B$,X$,Y$,Z$,N,M0,F$,A0,F0,S0,P,C 
WRITE  @D:X1,Y1,X,Y,M,L0,C0,A,B 
PRINT  @D, 2 : 

CLOSE 

REM 

PAGE 

GO  TO  5290 

ERROR  MESSAGE  IF  STORAGE  ATTEMPTED  BEFORE  PLOT  IS  READY 


REM 

PAGE 

PRINT  @32,21:30,50 

PRINT  "PLOT  NOT  COMPLETE — NOT  STOREABLE" 

HOME 

REM 

GO  TO  5540 
REM 

DELETE  Xl ,Y1 ,X,Y,M,L0 ,C0 ,A,B,R 

BRANCH  ON  DEVICE  TYPE 

GO  TO  1* (Dl=l )+2* (Dl=2  OR  Dl=3)  OF  5340,5430 

1.  RETRIEVAL  FROM  4907  FILE  MANAGER 


REM 

PRINT  "FILE  NAME:"; 

INPUT  W$ 

OPEN  W$ ; 1 , "R"  ,U$ 

READ  #1:N1,N2,N3,N4,A$,B$,X$,Y$,Z$,N,M0,F$,A0,F0,S0,P,C 

DIM  X1(N1) ,Y1(N3) ,X(N,M0) ,Y(N,M0) ,M(N) ,L0(N) ,C0(N) ,A(N) ,B(N) ,PJ 

READ  #1:X1,Y1,X,Y,M,L0,C0,A,B 

CLOSE  1 

GO  TO  5510 


2.  RETRIEVAL  FROM  4924  TAPE  DRIVE  OR  INTERNAL 
DRIVE  (NOTE:  4924  ASSUMED  AT  ADDRESS  2  OF 
GPIB;  SEE  LINE  9350) 


REM 

PRINT  "FILE  NUMBER:"; 
INPUT  D2 

D=D0*(D1=2)+33*(D1=3) 
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5470 

5480 

5490 

5500 

5510 

5520 

5530 

5540 

5550 


5560 

5570 

5580 

5590 

5600 

5610 

5620 

5630 

5640 

5650 

5660 

5670 

5680 

5690 

5700 

5710 

5720 

5730 


5740 

5750 

5760 

5770 

5780 

5790 

5800 


5810 

5820 

5830 

5840 

5850 

5860 

5870 

5880 

5890 

5900 


FIND  @D:D2 

READ  @D:Nl,N2,N3,N4,A$,B$,X$,Y$,Z$,N,M0,F$,A0,F0,S0,P,C 
DIM  Xl(Nl) ,Y1(N3) ,X(N,M0) ,Y(N,M0) ,M(N) ,L0(N) ,C0(N) ,A(N) ,B(N) ,R(N) 
READ  @D:X1,Y1,X,Y,M,L0,C0,A,B 
REM 

GOSUB  4410 
K=1 
REM 
RETURN 


USER  KEY  5:  CURVE  FITS 


PAGE 

IF  K ( 4 ) =0  THEN  6330 

PRINT  "CURVE  FITS  ( 0=NO , 1 =YES ) : " ; 

INPUT  C 

IF  C=0  THEN  6280 
IF  K ( 3 )  THEN  5630 
GOSUB  4150 
REM 
PRINT 

PRINT  "CURVE  FITS  AVAILABLE:" 

PRINT 

PRINT  "  1 .UNCONSTRAINED" 

PRINT  "  2. FIX  SLOPE" 

PRINT  "  3. FIX  INTERCEPT" 

PRINT  "  4. FIX  BOTH" 

PRINT 

PRINT  "CURVE  FIT  DESIRED:"; 

INPUT  C 

TRANSFORM  DATA  VIA  FUNCTIONS  DEFINING  PAPER  TYPE 

FOR  1=1  TO  N 

DELETE  X0 , Y0 , Z0 

DIM  X0(M(I) ) ,Y0(M(I) ) ,Z0(M(I) ) 

FOR  J=1  TO  M( I ) 

X0( J)=FNA(X( I,J) ) 

Y0(J)=FNB(Y(I,J) ) 

NEXT  J 

BRANCH  ON  CURVE  FIT  SELECTED 

GO  TO  C  OF  5820,5890,5940,6020 
REM 

Z0=X0  *Y0 

A( I ) =M ( I ) *SUM( Z0)-SUM(X0 )*SUM( Y0) 

Z0=X0*X0 

A(I)=A(I)/(M(I)*SUM(Z0)-SUM(X0)“2) 

B(I)=(SUM(Y0)-A(I)*SUM(X0))/M(l) 

GO  TO  6050 
REM 

PRINT  "SLOPE  OF  CURVE  NUMBER  "  ;  I ;  "  :  "  ; 
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5910 

5920 

5930 

5940 

5950 

5960 

5970 

5980 

5990 

6000 

6010 

6020 

6030 

6040 

6050 

6060 

6070 

6080 

6090 

6100 

6110 

6120 

6130 

6140 

6150 

6160 

6170 

6180 

6190 

6200 

6210 

6220 

6230 

6240 

6250 

6260 

6270 

6280 

6290 

6300 

6310 

6320 

6330 

6340 

6350 


INPUT  A( I ) 

B(I)=(SUM( YO)-A(I) *SUM(X0) )/M(l) 

GO  TO  6050 
REM 

PRINT  "INTERCEPT  OF  CURVE  NUMBER  "  ;  I ;  "  :  "  ; 

INPUT  B ( I ) 

Z0=X0  *Y0 

A( I )=SUM( ZO )-B( I ) *SUM(X0 ) 

zo=xo*xo 

A( I ) =A( I )/SUM( ZO ) 

GO  TO  6050 
REM 

PRINT  "SLOPE,  INTERCEPT  OF  CURVE  NUMBER  ";I;":"; 
INPUT  A ( I ) , B ( I ) 

REM 


CALCULATE  R-SQUARE 

Z0=SUM( Y0 )/M( I) 

Z0  =  Y0-  Z0 
Z0=Z0*Z0 
R(  I )  =  1/SUM ( ZO ) 

Z0=A( I ) *X0 
Z0=B( I )+Z0 
Z0=Y0-Z0 

zo=zo*zo 

R(  I ) =100* ( 1-SUM( ZO ) *R(  I )  ) 

NEXT  I 

PRINT  OUT  RESULTS  OF  CURVE  FITS 


PAGE 

FOR  1  =  1  TO  ( 3  4- ( N+3 )  ) /2 
PRINT 
NEXT  I 

PRI  USI  6210:"  CURVE  # " , "  SLOPE","  INTERCEPT","  R-SQUARE" 

IMAGE  8X , 4 ( 1 0 A, 5X ) 

PRINT 

FOR  1=1  TO  N 

PRINT  USING  "8X,5X,2D,3X,3(5X,7D.2D)":I,A(I) ,B(I) ,R(I) 

NEXT  I 
HOME 

GO  TO  6300 
REM 
PAGE 
REM 

K ( 5 ) =1 
GO  TO  6370 

ERROR  MESSAGE  IF  FIT  ATTEMPTED  BEFORE  DATA  INPUT 

REM 

PRINT  @32,21:25,50 

PRINT  "NO  DATA  INPUT  YET — CURVE  FIT  IMPOSSIBLE" 
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6360  HOME 
6370  REM 
6380  RETURN 

USER  KEYS  6  AND  7:  PLOT  ON  4051  OR  4662 


6390  PAGE  @D: 

CHECK  FOR  MISSING  PLOT  PARAMETERS 


6400 

IF  K  ( 1  ) 

THEN 

6420 

6410 

GOSUB 

100 

6420 

REM 

6430 

IF  K ( 2 ) 

THEN 

6450 

6440 

GOSUB 

1310 

6450 

REM 

6460 

IF  K (  3  ) 

THEN 

6480 

6470 

GOSUB 

4150 

6480 

REM 

6490 

IF  K  (  4  ) 

THEN 

6510 

6500 

GOSUB 

1410 

6510 

REM 

6520 

IF  K  (  5  ) 

THEN 

6540 

6530 

GOSUB 

5560 

6540 

REM 

SELECT  PEN  COLOR  AND  SPEED 

6550  PRINT  @D, 8 : F0 
6560  PRINT  @D, 32 : "BY" ; SO 

DEFINE  WINDOW  IN  UDU ' S 

6570  Ul=FNA(Xl ( 1 ) ) 

6580  U2=FNA( XI ( N1 ) ) 

6590  U3=FNB(Y1 ( 1 ) ) 

6600  U4=FNB ( Y1 ( N3 ) ) 

6610  WINDOW  Ul  ,U2  ,U3 ,U4 

DEFINE  VIEWPORT  IN  GDU ' S 

D 1  =  HORIZONTAL  LOCATION  OF  LOWER-LEFT  CORNER 
D2  =  HORIZONTAL  EXTENT 

D3  =  VERTICAL  LOCATION  OF  LOWER-LEFT  CORNER 
D4  =  VERTICAL  EXTENT 

6620  Dl=26* (D<>32 )+40* ( D=3 2 ) 

6630  D2  =  45*(DO32)+60*(D=32) 

6640  D3=l 5* ( D< >32 ) +20* ( D=32 ) 

6650  D4  =  45*(DO32)+60*(D=32) 

6660  VIEWPORT  Dl ,Dl+D2 ,D3 ,D3+D4 
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CALCULATE  SCALES  AND  TICK  LENGTHS 


6670  X2=(U2-U1)/D2 
6680  Y2=(U4-U3)/D4 
6690  X3=X2*(D2  MAX  D4)/75 
6700  Y3=Y2*(D2  MAX  D4)/75 

DEFINE  DEFAULT  HARDWARE  CHARACTER  SIZES 

6710  Cl=l .792*(D<>32 )+l .79* (D=32) 

6720  C2=l . 195* (D<>32)+1 . 195* (D=32) 

6730  C3  =  2 . 816 * ( D< >3 2 ) +2 . 6* ( D=3 2 ) 

6740  C4=l . 721 * ( D<>3 2 )  +  2 . 2* (D=3 2 ) 

DEFINE  RELATIVE  CHARACTER  SIZES  FOR  TICK  MARK  LABELS, 
AXIS  LABELS,  PLOT  TILE,  AND  LEGENDS 

6750  S1=0 .8* (D<>32 )+l* (D=32) 

6760  S2=1*(D<>32)+1* (D=32 ) 

6770  S3=0 .7* (D<>32 )+l* (D=32 ) 

DEFINE  OFFSETS  FOR  AXIS  LABELS  AND  PLOT  TITLE 

6780  Fl=7 * ( D<>3 2 ) +10 * ( D=32 ) 

6790  F2=10 *  ( DO 3 2  )  +  14 *  ( D=32  ) 

6800  F3=7* (D<>32 )+9* (D=32 ) 

BRANCH  ON  FRAME  TYPE  (ZERO  MEANS  NO  FRAME  OR  LABELS) 

6810  IF  A0  =  0  THEN  8060 

6820  GO  TO  A0  OF  6830,6990,7190 

1.  DRAW  AXES 

6830  REM 

6840  MOVE  @D:U1,U4 

6850  FOR  I=N3  TO  2  STEP  -1 

6860  D2=FNB ( Yl ( I ) ) 

6870  DRAW  @D:U1,D2 

6880  DRAW  @D:Ul+X3,D2 

6890  DRAW  @D:Ul,D2 

6900  NEXT  I 

6910  DRAW  @D:U1 ,U3 

6920  FOR  1=2  TO  N1 

6930  D1=FNA(X1 ( I ) ) 

6940  DRAW  @D:D1,U3 

6950  DRAW  @D:D1,U3+Y3 

6960  DRAW  @D:Dl,U3 

6970  NEXT  I 

6980  GO  TO  7570 
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6990 

7000 

7010 

7020 

7030 

7040 

7050 

7060 

7070 

7080 

7090 

7100 

7110 

7120 

7130 

7140 

7150 

7160 

7170 

7180 

7190 

7200 

7210 

7220 

7230 

7240 

7250 

7260 

7270 

7280 

7290 

7300 

7310 

7320 

7330 

7340 

7350 

7360 

7370 

7380 

7390 

7400 

7410 

7420 

7430 

7440 

7450 

7460 

7470 


2.  DRAW  GRID 


REM 

MOVE  @D:U1,U3 
DRAW  @D:U1,U4 
D2=U4 

FOR  1=2  TO  Nl 
Dl=FNA(Xl ( I ) ) 

DRAW  @D:D1,D2 
D2=U3* (D2=U4 )+U4* (D2=U3 ) 
DRAW  @D:D1,D2 
NEXT  I 

DRAW  @D:U2,U3 
DRAW  @D:U1,U3 
Dl  =  Ul 

FOR  1=2  TO  N3 
D2=FNB(Y1 ( I ) ) 

DRAW  @D:D1,D2 
Dl=Ul * ( Dl=U2 ) +U2* ( Dl=Ul ) 
DRAW  @D:D1,D2 
NEXT  I 
GO  TO  7570 

3.  DRAW  BORDER 


REM 

MOVE  @D:U1,U3 
IF  Nl  =  2  THEN  7280 
FOR  1=2  TO  Nl-1 
Dl=FNA( Xl ( I ) ) 
DRAW  @D:D1,U3 
DRAW  @D: Dl , U3+Y3 
DRAW  @D:D1,U3 
NEXT  I 
REM 

DRAW  @D:U2,U3 
IF  N3=2  THEN  7370 
FOR  1=2  TO  N3-1 
D2=FNB(Y1 (I ) ) 
DRAW  @D:U2,D2 
DRAW  @D: U2-X3 , D2 
DRAW  @D:U2,D2 
NEXT  I 
REM 

DRAW  @D:U2,U4 
IF  Nl=2  THEN  7460 


FOR  I=N1— 1  TO  2  STEP  - 
D1=FNA(X1 ( I ) ) 

DRAW  @D : Dl , U4 
DRAW  @D:D1,U4-Y3 
DRAW  @D:D1,U4 
NEXT  I 
REM 

DRAW  @D:Ul,U4 


1 
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7480 

7490 

7500 

7510 

7520 

7530 

7540 

7550 

7560 

7570 

7580 

7590 

7600 

7610 

7620 

7630 

7640 

7650 

7660 

7670 

7680 

7690 

7700 

7710 

7720 

7730 

7740 

7750 

7760 

7770 

7780 

7790 

7800 

7810 

7820 

7830 

7840 

7850 

7860 

7870 

7880 

7890 


IF  N3=  2  THEN  7550 

FOR  I=N3-1  TO  2  STEP  -1 
D2=FNB( Yl ( I ) ) 

DRAW  @D:Ul,D2 
DRAW  @D : U1+X3 , D2 
DRAW  @D:U1,D2 
NEXT  I 
REM 

DRAW  @D: U1 ,U3 
REM 


PRINT  X-AXIS  LABEL 

PRINT  @D,17:S1*C1,S1*C3 
PRINT  @D , 25 : 0 

MOVE  @D:  ( U1+U2-S1 *X2 * (C1*(LEN(X$)-1)+C2) )/2,U3-Fl*Y2 
PRINT  @D:X$ 

PRINT  Y-AXIS  LABEL 

IF  D=32  THEN  7670 
PRINT  @D, 25:90 

MOVE  @D:U1-(F2-S1*C4 ) *X2 , ( U3+U4-S1 * Y2* ( Cl  * ( LEN ( Y$ ) -1 )+C2) )/2 
PRINT  @D:Y$ 

GO  TO  7730 
REM 

FOR  1=1  TO  LEN ( Y$ ) 

MOVE  Ul-F2*X2 ,  ( U3+U4  +  Y2 *Sl * ( C3* ( LEN ( Y$ ) -1 )+C4 ) )/2-Y2* ( C4+C3* ( 1-1 ) ) 
W$=SEG( Y$ ,1,1 ) 

PRINT  W$ 

NEXT  I 
REM 


PRINT  AND  UNDERLINE  PLOT  TITLE 

IF  LEN ( Z$ )=0  THEN  7830 
PRINT  @D, 17 : S2*C1 ,S2*C3 
PRINT  @D, 25 : 0 

MOVE  @D:  ( Ul+U2-S2*X2 * (C1*(LEN( Z$ )  — 1 )+C2) )/2,U4+F3*Y2 
PRINT  @D:Z$ 

MOVE  @D , 21 : ( Ul+U2+S2*X2 * ( Cl* ( LEN ( Z$ )-l )+C2 ) )/2 ,U4+( F3-0 . 2*S2*C3 ) *Y2 
MOVE  @D,20: ( Ul+U2-S2*X2 * ( Cl  * ( LEN( Z$ ) -1 ) +C2 ) )/2 ,U4  +  ( F3-0 . 2*S2*C3 ) *Y2 
MOVE  @D,21:  ( U1+U2-S2 *X2 * ( C 1 * ( LEN ( Z$ ) - 1 ) +C2 )  ) /2 , U4+ ( F3-0 . 4 *S2*C3 ) *Y2 
MOVE  @D, 20 : ( U1  +  U2+S2 *X2* ( Cl  * ( LEN ( Z$ ) -1 )+C2)  )/2 , U4  +  ( F3-0 . 4 *S2*C3 ) *Y2 
REM 

PRINT  @D,17:S3*C1,S3*C3 
PRINT  @D, 2  5 : 0 

PRINT  TICK  MARK  LABELS  ON  X  AXIS 

IF  N2=0  THEN  7950 
U$=A$ 

FOR  1=1  TO  N2 

W$=SEG ( U$ , 1 , POS (U$,"#M ,1)-1) 
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7900 

7910 

7920 

7930 

7940 

7950 


7960 

7970 

7980 

7990 

8000 

8010 

8020 

8030 

8040 

8050 

8060 

8070 

8080 

8090 


8100 


8110 


8120 


8130 

8140 


8150 

8160 

8170 

8180 

8190 

8200 

8210 

8220 


U$=SEG ( U$ , POS ( U$ , " # " , 1 )+l , LEN ( U$ ) -POS ( U$ , ” # " , 1 ) ) 

Dl=FNA ( VAL ( W$ ) ) 

MOVE  @D: Dl-S3*X2 * ( Cl* ( LEN ( W$ ) -1 ) +C2 )/2 , U3-S3*Y2* ( 1 . 25*C4 ) 
PRINT  @D:W$ 

NEXT  I 
REM 


PRINT  TICK  MARK  LABELS  ON  Y  AXIS 

IF  N4=0  THEN  8050 
U$  =  B$ 

FOR  1=1  TO  N4 

W$  =  SEG  ( U$  ,  1 ,  POS  (  U$  ,  "  # 11  ,1)-1) 

U$=SEG(U$,POS(U$, "#" ,1 )+l ,LEN(U$ )-POS(U$, "#" ,1 ) ) 

Dl=FNB ( VAL( W$ ) ) 

MOVE  @D:U1-S3*X2*(C1*LEN(W$)-0.2*(C1-C2) ) , Dl-S3*Y2* ( 1 . 125 *C4 ) /2 
PRINT  @D:W$ 

NEXT  I 
REM 
REM 

PRINT  @D,17:S3*C1 ,S3*C3 

PRINT  @D, 25 : 0 

U$=F$ 


LOOP  OVER  N  DATA  CURVES 
FOR  1=1  TO  N 

1.  SELECT  COLOR  FOR  I ' TH  CURVE 
PRINT  @D, 8 : CO ( I ) 

2.  BRANCH  ON  CURVE  TYPE 

GO  TO  1*(L0(I)<=3)+2*(L0(I)=4)+3*(L0(I)=>5)  OF  8130,8710,8810 
2. A  DASHED  LINES  (CURVE  TYPE  1  THROUGH  3) 


REM 

IF  M( I ) <= 1  THEN  8690 

2.A.1  LOAD  IN  DASHED  LINE  PARAMETERS  (IN  GDU'S) 
(NOTE:  1 1  ST  DASHED  LINE  HAS  3  ELEMENTS, 

2 ' ND  HAS  5,3'RD  HAS  3) 


RESTORE  8160 

DATA  3, 0.2, 0.6,0. 2,5,0.2,0.6,0.01,0.6,0.2,3,0.005,0.6,0.005 
FOR  J=1  TO  L0 ( I ) 

READ  L 
DELETE  P0 
DIM  P0(L) 

READ  P0 
NEXT  J 
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2. A. 2  TRANSFORM  DATA 


8230 

8240 

8250 

8260 

8270 

8280 

DELETE  X0,Y0 

DIM  X0 ( M ( I ) ) ,  YO  (  M(  I ) ) 

FOR  J=1  TO  M( I ) 

XO ( J)=FNA(X(I , J) ) 

YO ( J ) =FNB( Y ( I, J) ) 

NEXT  J 

2. A. 3  CALCULATE  CURVE  LENGTH  IN  GDU'S 

8290 

8300 

*  8310 

8320 

T0  =  0 

FOR  J=2  TO  M( I ) 

T0=T0+SQR( ( (XO(J)-XO(J-l) )/X2)-2+( ( YO ( J )-Y0 ( J-l ) )/Y2)=2) 
NEXT  J 

2. A. 4  CALCULATE  ADJUSTED  DASHED  LINE  PARAMETERS 

8330 

8340 

Ml=l  MAX  I NT (T0/SUM(P0)+0.5) 

P0=T0/(M1*SUM(P0 ) )*P0 

2. A. 5  MOVE  TO  START  OF  CURVE 

8350 

8360 

8370 

8380 

J=1 

D1=X0 ( J ) 

D2  =  YO ( J  ) 

MOVE  @D:Dl,D2 

2. A. 6  LOOP  OVER  Ml  PATTERNS  AND  L  ELEMENTS  IN 
EACH  PATTERN 

8390 

8400 

8410 

8420 

8430 

FOR  11=1  TO  Ml 

FOR  Jl=l  TO  L 

D3  =  0 

D4=0 

REM 

<3 

8440 

8450 

8460 

8470 

8480 

8490 

8500 

8510 

8520 

8530 

8540 

2. A. 7  CALCULATE  DISTANCE  TO  NEXT  DATA  POINT; 

IF  GREATER  THAN  PATTERN  ELEMENT,  INTER¬ 
POLATE  TO  FIND  WHERE  TO  GO  AND  SET  FLAC- 
TO  INDICATE  DONE  WITH  ELEMENT — IF  LESS, 

GO  TO  NEXT  DATA  POINT 

D5= (X0 ( J+l )-Dl )/X2 

D6=(Y0(J+1)-D2)/Y2 

D7=SQR( D5=2+D6=2 ) 

IF  D4+D7=>P0 ( Jl )  THEN  8530 

J=J  +  1 

D1=X0 ( J  ) 

D2  =  Y0 ( J  ) 

D4  =  D4  +  D7 

GO  TO  8590 

REM 

IF  D7<=1 . OE-4  THEN  8570 
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8550 

8560 

8570 

8580 

8590 

8600 

8610 

8620 

8630 

8640 

8650 

8660 

8670 

8680 

8690 

8700 

8710 

8720 

8730 

8740 

8750 

8760 

8770 

8780 

8790 

8800 

8810 

8820 

8830 

8840 

8850 

8860 

8870 

8880 

8890 

8900 

8910 


D1=D1+X2*(P0(J1)-D4)*D5/D7 
D2=D2+Y2* ( P0 ( J 1 ) -D4 ) *D6/D7 
REM 
D3  =  1 
REM 

2. A. 8  IF  EVEN-NUMBERED  ELEMENT  BEING  PROCESSED, 

MOVE;  IF  ODD-NUMBERED,  DRAW 

IF  J 1  =  2 * INT ( Jl/2 )  THEN  8630 

DRAW  @D:D1,D2  * 

GO  TO  8650 
REM 

MOVE  @D:D1,D2 
REM 


2. A. 9  TEST  IF  DONE  WITH  CURRENT  ELEMENT 

IF  D3<>1  THEN  8430 
NEXT  J1 
NEXT  11 
REM 

GO  TO  9180 

2.B  SOLID  LINE  (CURVE  TYPE  4) 


REM 

Dl=FNA( X ( I , 1 ) ) 

D2=FNB(Y( 1,1 ) ) 

MOVE  @D:  Dl  , D2 
FOR  J=1  TO  M ( I ) 

Dl  =  FNA( X ( I , J)  ) 

D2=FNB( Y( I , J ) ) 

DRAW  @D: Dl , D2 
NEXT  J 
GO  TO  9180 

2.C  SYMBOLS  (CURVE  TYPES  5  THROUGH  9) 


REM 

FOR  J=1  TO  M ( I ) 

Dl=FNA(X ( I , J ) ) 

D2=FNB ( Y ( I , J ) ) 

GO  TO  L0 ( I ) -4  OF  8860,8920,8980,9050,9110 
2.C.1  CROSSES 


REM 

MOVE  @D: Dl+1 . 5*X3 ,D2 
DRAW  @D:D1-1.5*X3,D2 
MOVE  @D: Dl ,D2-1 . 5*Y3 
DRAW  @ D: Dl , D2  +  1 . 5 *Y3 
GO  TO  9160 
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2.C.2  X-SHAPED 


8920 

8930 

8940 

8950 

8960 

8970 


8980 

8990 

9000 

9010 

9020 

9030 

9040 


9050 

9060 

9070 

9080 

9090 

9100 


9110 

9120 

9130 

9140 

9150 

9160 

9170 

9180 


9190 

9200 

9210 

9220 

9230 

9240 


9250 

9260 

9270 


REM 

MOVE  @D: D1+X3 , D2+Y3 
DRAW  @D: D1-X3 ,D2-Y3 
MOVE  @D: D1-X3 , D2+Y3 
DRAW  @D: D1+X3 ,D2-Y3 
GO  TO  9160 

2.C.3  SQUARES 


REM 

MOVE  @D: D1-X3 , D2-Y3 
DRAW  @D: D1-X3 , D2+Y3 
DRAW  @D: D1+X3 , D2+Y3 
DRAW  @D: D1+X3 , D2-Y3 
DRAW  @D: D1-X3 , D2-Y3 
GO  TO  9160 

2.C.4  TRIANGLES 


REM 

MOVE  @D:D1 , D2+2* Y3/SQR( 3 ) 
DRAW  @D:Dl+X3,D2-Y3/SQR(3) 
DRAW  @D: D1-X3 , D2-Y3/SQR( 3 ) 
DRAW  @D: Dl ,D2+2*Y3/SQR(3 ) 
GO  TO  9160 


2.C.5  UPSIDE-DOWN  TRIANGLES 


REM 

MOVE  @D:D1 ,D2-2*Y3/SQR( 3 ) 
DRAW  @D:Dl-X3,D2+Y3/SQR(3) 
DRAW  @D: D1+X3 , D2+Y3/ SQR( 3 ) 
DRAW  @D: Dl , D2-2*Y3/SQR( 3 ) 
REM 

NEXT  J 
REM 


3.  PRINT  LEGEND  FOR  I ' TH  CURVE 
W$=SEG(U$,1 ,POS(U$, "#" ,1 )-l) 

U$=SEG(U$,POS(U$,"#" , l)+l,LEN(U$)-POS(U$,"#",l) ) 
Dl=FNA(X( I ,M( I ) ) ) 

D2=FNB ( Y ( I , M ( I ) ) ) 

MOVE  @D:D1+3*X3,D2-0.5*S3*Y2*(1.125*C4) 

PRINT  @D:W$ 


4.  PLOT  CURVE  FIT  IF  REQUIRED 

IF  C=0  THEN  9280 

MOVE  @D:U1 ,A( I ) *Ul+B( I ) 

DRAW  @D:U2,A(I)*U2+B(I) 
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9280  REM 
9290  NEXT  I 

RETURN  PEN  TO  HOLDER  AND  EXIT 

9300  PRINT  @D,8:0 
9310  HOME  @D: 

9320  RETURN 

INITIALIZATION  ROUTINE 

ARRAY  K  IS  USED  TO  KEEP  TRACK  OF  WHICH  USER  KEYS 
HAVE  BEEN  DEPRESSED;  LINE  9350  DEFINES  THE  ADDRESS  OF 
THE  4924  TAPE  DRIVE 

9330  DIM  K(5 ) ,R$ (1) 

9340  R$=REP("",1,0) 

9350  D0=2 

9360  IF  LEN ( R$ ) > 0  THEN  9390 
9370  K=0 
9380  R$= "  " 

9390  REM 
9400  RETURN 
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